UNCLASSIFIED 


Reproduced 
luf  ike 

ARMED  SERVICES  TECHNICAL  INFORMAnON  AGENCT 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Gkjvemment  thereby  Incurs  no  responsibility,  nor  any 
obligation  \diatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  coiT>oration,  or  conveying  any  rl^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


KVpott  190/ 


v'iH 

wUhk&.'i 


1 

tH  «ktpW>ff.>^.4UU0W— 1 

^KiTr 

Wi  m 

iiLi* 


COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR 
MARINE  CONTROL-SURFACE  FLUTTER 


mmm& 


RolHt  C.  Loibowiti  «mI  DmoIcI  J.  Bolt 


smciUttL 


APPLIED 

MADflUTICS 


ASTI  A 

EiSEXH  nr?.. 

KB4  1963 
mk 


STRUCTURAL  MECHANICS  LABORATORY 
RESEARCH  AND  DEVELOPMENT  REPORT 


August  1962 


Roport  1567 


COMPARISON  OF  THEORY  AND  EXPERIMENT  FOR 
MARINE  CONTROL-SURFACE  FLUTTER 


Ralph  C.  Laibowitz  and  Donald  J.  Btiz 


Rtprinted  from  a  pop«r  prasentad  at  the  Fourth  Symposium  on 
In'ovoI  Hydrodynamics  on  "Ship  Propulsion  and  Hydroelosticity" 
at  Washington,  D.C.,  27—31  August  1962,  Vol.  3  of  ACR-73, 
Office  of  Naval  Research.  Discussion  is  omitted.  Symposium 
sponsored  by  ONR,  Department  of  the  Navy,  and  ossisted  by 
Webb  Institute  of  Novai  Architecture. 


August  1962 


Report  1567 


Leibowitz  and  Belz 


TABLE  OF  CONTENTS 

Page 

Abstract  . 

Introduction  .  713 

Background .  714 

TMB  Control  Surface  Flutter  Apparatus  .  715 

Flutter  Analysis .  717 

Equations  of  Motion .  717 

Digital  Computer  Determination  of  Damping  Ratios  .  718 

Analog  Determination  of  Flutter  Response .  721 

Experimental  Work .  725 

Analytical  and  Experimental  Results  .  725 

Discussion  of  Results .  727 

Conclusions .  746 

Recommendations .  747 

Acknowledgments  .  749 

Appendix  A  Comparison  of  Predicted  and  Observed 

Damping  Ratios  and  Associated  Frequen¬ 
cies  for  the  TMB  Control  Surface  Flutter 

Apparatus  .  750 

Appendix  B  Numerical  Values  of  the  Parameters  mh,  L, 

S,  and  .  820 

Appendix  C  Relationship  of  the  Modal  Frequencies  to  the 

Phase  Angle  between  Y  and  9 . *  •  •  821 

\ppendix  D  Relationship  between  the  Parameter  L  and 

the  Analytical  Expressions  for  Hydrodynamic 

Lift  and  Moment .  824 

Appendix  E  Cr  teria  for  Selection  of  Initial  Conditions  which 
Allow  One  Mode  to  Predominate  in  the  Total 

Response .  826 

Appendix  F  The  Effect  of  Initial  Conditions  on  the  Vibration 

of  a  Ship's  Rudder .  830 

Appendix  G  Analogy  between  TMB  Control  Surface  Flutter 
Apparatus  and  the  Rudder  Configuration  of 

USS  FORREST  SHERMAN  (DD  931) .  831 

Appendix  H  The  Classical  Theodorsen  Equations  and  their 

Questionable  Applicability  to  Rudders .  833 

References  .  846 

Bibliography  .  847 


706 


Leibowitz  and  Belz 


LIST  OF  FIGURES 


Page 

Figure  1 

Schematic  Representation  of  the  TMB  Control 
Surface  Flutter  Apparatus 

716 

Figure  2 

Sample  Format  of  Digital  Computer  Output 

Data 

722 

Figure  3 

Analog  Response  Showing  Stable  and  Unstable 
Oscillations 

724 

Figure  4 

Comparison  of  Computed  and  Experimental 
Damping  Ratios  (Configuration  TMB  1222, 
mh  =  0.  0  lb-sec2) 

728 

Figure  5 

Comparison  of  Computed  and  Experimental 
Damping  Ratios  (Configuration  TMB  1222, 
mh  =  0.  3  lb-sec^) 

729 

Figure  6 

Comparison  of  Computed  and  Experimental 
Damping  Ratios  (Configuration  TMB  1222,  mh  = 
2.  0  lb-sec^) 

730 

Figure  7 

Comparison  of  Computed  and  Experimental 
Damping  Ratios  (Configuration  B,  mh  =  3.  53 
lb-sec2) 

731 

Figure  8 

Comparison  of  Computed  and  Experimental 
Damping  Ratios  (Configuration  B,  mh  =  3.  68 
lb-sec^) 

732 

Figure  9 

Comparison  of  Computed  and  Experimental 
Damping  Ratios  (Configuration  B,  mh  =  3.  82 
lb-sec2) 

733 

Figure  10 

Comparison  of  Computed  and  Experimental 
Damping  Ratios  (Configuration  B,  mh  =  3.  94 
lb-sec^) 

734 

Figure  1 1 

Comparison  of  Computed  and  Experimental 
Damping  Ratios  (Configuration  B,  mh  =  4,  05 
lb-sec^) 

735 

Figure  12 

Comparison  of  Computed  and  Experimental 
Damping  Ratios  (Configuration  B,  mh  =  4.  21 
lb-sec^) 

736 

Figure  13 

Comparison  of  Computed  and  Experimental 
Damping  Ratios  (Configuration  B,  mh  =  4.  32 
lb-sec2) 

737 

Figure  14 

Comparison  of  Computed  and  Experimental 
Damping  Ratios  (Configuration  D,  mh  =  4,  05 
lb-sec2) 

738 

Figure  15 

Comparison  of  Computed  and  Experimental 
Damping  Ratios  (Configuration  D,  mh  =  4.  21 
lb-sec2) 

739 

707 


Leibowitz  and  Belz 


Figure  16  Comparison  of  Computed  and  Experimental  Damp¬ 
ing  Ratios  (Configuration  D,  mh  =  4.  32  lb-sec^)  740 

Figure  17  Comparison  of  Computed  and  Experimental  Criti¬ 
cal  Flutter  Speeds  for  an  Assumed  Value  of  L 
(Configuration  B) 

Figure  18  Comparison  of  Computed  and  Experimental  Criti¬ 
cal  Flutter  Speeds  for  an  Assumed  Value  of  L 
(Configuration  D)  742 

APPENDIX  A 

Figures  19  through  86  -  Comparison  of  Predicted  and  Ob¬ 
served  Damping  Ratios  and  Associated  Frequencies  for  the 
TMB  Control  Surface  Flutter  Apparatus  750 

APPENDIX  C 

Figure  87  Sign  Convention  for  Two-Degree-of-Freedom- 

Mass  Spring  System  and  Hydrofoil  System  822 

APPENDIX  H 

Figure  88  Forces  and  Moment  on  Foil  Moving  in  a  Uniform 

Stream  834 


708 


Leibowitz  and  Belz 


LIST  OF  TABLES 


Page 


Table  1  Comparison  of  Notation  Used  in  Figure  2  with 

Notation  of  this  Report  .  723 

Table  2  Parameters  of  the  TMB  Control  Surface  Flutter 

Apparatus  for  Various  Configurations  .  .  . .  726 

Table  3  Numerical  Values  of  Parameters  mh,  L,  S,  and 
no  1  for  which  Analytical  Frequency  and  Stabil¬ 
ity  Predictions  Have  Been  Made .  751 

Table  4  Relation  between  Frequency  and  Phase  Angle 
for  Common  Two- Mass  System  and 
Hydrofoil  of  Flutter  Apparatus .  823 


709 


Cd  W  to 


Leibowitz  and  Belz 


Symbol 
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al 


10*  ®11 
12.  B21 

22 
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C 


C 


c 


C{k) 

c 


e 


h 


I 


j 


NOTATION 

Definition  Dimensions 

Lift  constant  of  the  hydrofoil  (lift  force  lb-sec^ 
per  unit  attack  angle  per  unit  velocity  in.  ^ 
squared) 

Distance  from  midchord  to  axis  of  rota-  in. 
tion  (positive  aft  of  midchord) 

Symbols  for  coefficients  of  the  hydrofoil  varying 
equations  of  motion 


Semichord  length  of  hydrofoil  in. 

Linearized  damping  constant  for  the  lb- sec /in, 

translational  degree  of  freedom  at  zero 

speed 

Critical  damping  constant  (general  sym-  in.  -lb-sec(rot.  ) 
bol  employed  for  rotational  or  transla-  lb- sec /in. 

tional  degree  of  freedom)  (trans.  ) 

Theodorsen  function  Dimensionless 


Linearized  damping  constant  for  the  in.  -lb-sec 

rqtaUonal  degree  of  freedom  at  zero 

speed 

Base  of  Natural  Logarithms  (2.  718.  .  .  )  Dimensionless 


Oscillatory  lift  force  less  the  added  mass  lb 
and  zero  hydrodynamic  damping  effects 

Disteuice  from  the  axis  of  rotation  to  the  in. 
c.  g.  of  the  rotating  assembly,  based  on 
mass  plus  added  mass  (positive  if  c.  g. 
is  downstream) 

Effective  mass  moment  of  inertia  of  the  in.  -lb- sec 
rotating  assembly  with  respect  to  its 
axis,  including  the  added  mass  moment 
of  inertia 

The  square  root  of  minus  one  (  V-1  )  Dimensionless 
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Symbol 

K 

Definition 

Translational  spring  constant 

Dimensions 
lb /in. 

k 

Torsional  spring  constant 

in.  -lb 

L 

Distance  of  the  center  of  lift  from 
the  axis  of  rotation  (positive  if  center 
of  lift  is  forward  of  the  axis) 

in. 

M 

Mass  of  the  apparatus  that  moves  only 
in  translation 

Ib-sec^/in. 

Hydrodynamic  oscillatory  moment  ex¬ 
erted  on  the  hydrofoil  about  the  axis  of 
rotation  less  the  added  mass  and  zero 
speed  hydrodynamic  damping  effects 

in.  -lb 

m 

Mass  of  the  hydrofoil  including  the  mass 
of  the  entire  assembly  that  rotates  with 
it  and  the  added  mass  for  translation 

Ib-sec^/in. 

m 

Added  mass  of  the  hydrofoil  in  trans¬ 
lation  or  Y  degree  of  freedom 

Ib-sec2/in. 

n 

Multiplier  of  t  in  assumed  solution 
d  m  e*'^  (See  p.  827  in  Appendix  E) 

sec"  ^ 

^10’ Jll 
Ri2»^20 
R21>  ^22 

Symbols  for  coefficients  of  the  hydro¬ 
foil  equations  of  motion. 

varying 

S 

Hydrofoil  towing  speed 

in.  /sec  in  com¬ 
putations,  con¬ 
verted  to  knots 
for  data  presen¬ 
tation 

Sa 

Moment  about  the  forward  quarter- 
point  of  the  rotational  structural  mass 

lb-sec^ 

T 

Period  of  an  oscillation 

sec 

t 

Elapsed  time 

sec 

*0 

Duration  of  an  impulsive  force 

sec 
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Symbol  Definition  Dimensions 

tj,  Successive  numerical  values  of  sec 

elapsed  time 


Y  Translational  displacement  of  the  hydro-  in. 

foil  axis  normal  to  the  flow  and  to  the 
axis  of  rotation  (measured  from  the 
equilibrium  position) 

Yq  Numerical  value  of  Y  at  t  =  0  sec  in, 

oq  Preset  angle  of  attack  of  hydrofoil  degrees 

5  Logarithmic  decrement  dimensionless 


6  Rotational  displacement  of  the  control  radians 

surface  from  the  direction  of  flow 


X 

X* 

y- 


CJ 


X 


a 


Numerical  value  0  at  t  =  0  sec 

radians 

Complex  circular  frequency  of  vibration 

sec"  ^ 

Complex  conjugate  of  X 

sec"  ^ 

Real  part  of  X  indicating  the  degree  of 
damping 

sec"  ^ 

Ratio  of  circumference  to  diameter  of 
a  circle  (3.  141. ,  .  ) 

dimensionless 

Phase  angle  separating  Y  and  9  re¬ 
sponses 

degrees 

Circular  frequency  of  vibration  (mag¬ 
nitude  of  imaginary  part  of  X  ) 

sec"  ^ 

The  distance  of  the  center  of  mass 
aft  of  the  rotational  axis 

dimensionless 
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COMPARISON  OF  THEORY  AND  EXPERIMENT 

FOR 

MARINE  CONTROL- SURFACE  FLUTTER 


Ralph  C.  Leibowitz  and  Donald  J.  Belz 
David  Taylor  Model  Basin 
Structural  Mechanics  Laboratory 

ABSTRACT 

Both  the  Extended  Simplified  Flutter  Analysis  and  Modified 
Theodorsen  Flutter  Analysis,  proposed  by  McGoldrick  and  Jewell, 
are  applied  to  the  TMB  CONTROL  SURFACE  FLUTTER  APPARA¬ 
TUS.  Predictions  of  vibrational  stability  and  instability  based  on 
these  analyses  are  compared  with  stable  and  unstable  (classical 
flutter)  vibrations  observed  in  the  apparatus  for  towing  speeds  in  the 
range  of  0  to  20  knots. 

The  comparison  of  predicted  and  observed  values  of  (1) 
damping  and  (2)  critical  flutter  speeds  shows  that  the  Modified  Theo¬ 
dorsen  Analysis  gives  a  consistently  better  agreement  with  experi¬ 
mental  data  than  does  the  Extended  Simplified  Analysis;  moreover, 
the  results  of  the  former  analysis  are  in  good  agreement  with  avail¬ 
able  experimental  data,  whereas  the  results  of  the  latter  analysis 
are  not. 


To  extend  the  range  of  mass  unbalance,  speed,  and  other 
parameters  that  show  good  agreement  between  theory  and  experi¬ 
ment,  certain  studies  that  will  yield  refinements  to  the  Modified 
Theodorsen  Analysis  are  proposed. 


INTRODUCTION 

The  possibility  of  flutter  in  marine  control  surfaces  came 
pointedly  to  the  attention  of  naval  research  personnel  when  sea 
trials  of  USS  FORREST  SHERMAN  (DD  931)  revealed  that  severe 
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vibrations  were  transmitted  to  the  hull  by  the  rudders.  ^  '  This  indi¬ 
cated  the  need  for  a  thorough  exploration  of  the  possibility  of  marine 
control- surface  flutter;  for  it  may  be  possible  that  as  the  speeds  of 
other  ships  are  increased  flutter  may  occur  more  commonly  within 
the  range  of  operating  speeds.  For  this  reason,  methods  of  pre¬ 
dicting  flutter  warrant  serious  consideration  in  the  field  of  naval 
architecture. 

The  present  comparison  of  theoretical  and  experimental 
results  is  intended  as  a  help  to  the  researchers  under  contract  to 
the  Bureau  of  Ships  in  their  search  for  a  suitable  method  of  flutter 
prediction.  Such  a  method  will  permit  ship  designers  to  establish 
stability  margins  of  safety  from  flutter  for  ships  potentially  capable 
of  control- surface  flutter 

In  this  report,  the  stability  of  a  TMB  Control  Surface 
Flutter  Apparatus  has  been  analyzed  by  computing  on  a  digital  com¬ 
puter  the  damping  associated  with  the  apparatus;  stability  predic¬ 
tions  based  on  the  se  computations  are  compared  with  available  ex¬ 
perimental  data.  These  data  include  the  results  of  tests  performed 
on  the  flutter  apparatus  published  in  Reference  2  and  recently  pub¬ 
lished  results  of  similar  tests  recently  performed  by  the  TMB 
Hydromechanics  Laboratory.  The  recent  data  are  of  special  in¬ 
terest  because,  as  a  result  of  chemges  made  in  certain  physical 
parameters  of  the  apparatus,  classical  flutter***  was  observed  as 
a  hydroelastic  phenomenon. 

In  addition  to  the  digital  computations,  an  analog  solution 
of  the  equations  of  motion^  used  in  the  analysis  of  the  flutter 
apparatus  was  undertaken.  Qualitative  information  on  the  total 
motional  response  was  obtained  from  this  solution. 

BACKGROUND 

A  Control  Surface  Flutter  Apparatus,  shown  schematically 
in  Fig.  1,  was  built  and  tested  at  the  David  Taylor  Model  Basin  as 


''^‘References  are  listed  on  page  846 

*'^A  preliminary  report  on  the  present  study  was  published  in  letter 
form  as  Reference  3. 

***Classical  flutter  is  a  dynamically  unstable,  self- excited  vibra¬ 
tion  of  an  oscillatory  system  immersed  in  a  field  of  fluid  flow.  ^ 
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a  joint  effort  of  the  Structural  Mechanics  and  Hydromechanics 
Laboratories.  ^  Experimental  measurements  of  overall  damping’'' 
associated  with  the  hydrofoil  of  the  apparatus  were  made  from  data 
obtained  while  the  apparatus  was  towed  with  its  foil  submerged. 

In  addition,  two  auialytic  approaches  to  the  problem  of  pre¬ 
dicting  flutter  were  devised  and  designated  as  the  Extended  Simpli¬ 
fied  Analysis  and  the  Modified  Theodor  sen  Analysis.  ^ 

The  characteristic  equation  obtained  by  expansion  of  the  determinant 
of  the  coefficients  of  the  differential  equations  of  motion  was  con¬ 
verted  to  algebraic  form  for  both  approaches,  and  the  Routh  dis¬ 
criminant  of  these  characteristic  equations  was  determined  in  order 
to  predict  the  stability  of  the  hydrofoil.  The  results  of  this  study 
were  published  in  Reference  2. 

That  reference  stated  that  analytical  predictions  of  in¬ 
stability,  based  on  the  Routh  criteria  alone,  were  not  verified  ex¬ 
perimentally.  However,  marked  reductions  in  the  damping  of 
stable  oscillations  were  observed.  This  associated  condition  of 
barely  stable  vibration  was  defined  as  subcritical  flutter.  Because 
subcritical  flutter  is  generally  not  predictable  from  the  Routh  cri¬ 
teria,  the  present  study  was  initiated  to  compute  damping  directly 
in  order  to  furnish  a  more  general  prediction  of  stability  for  the 
flutter  apparatus. 

TMB  CONTROL  SURFACE  FLUTTER  APPARATUS 

The  TMB  Control  Surface  Flutter  Apparatus**  was  designed 
to  provide  a  means  of  "checking  flutter  euialyses  based  on  various 
assumptions  as  to  the  nature  of  the  oscillatory  lift  forces  and 
moments.  "2 

The  mechanism  consists  of  a  relatively  rigid  NACA  0015- 
section  hydrofoil  so  suspended  that  it  rotates  about  a  vertical  axis 
and  translates  normal  to  the  path  of  a  towing  carriage  upon  which 
the  Control  Surface  Flutter  Apparatus  is  mounted.  Deflections  in 


^Overall  damping,  as  defined  on  page  4  of  Reference  2,  includes 
all  effects  contributing  to  the  rate  of  decay  or  buildup  of  a  free 
vibration. 

**A  detailed  description  of  this  apparatus  may  be  found  in  Reference 

2. 
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Figure  1  -  Schematic  Representation 
of  the  TMB  Control  Surface  Flutter 
Apparatus 
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the  two  degrees  of  freedonn,  0  (rotation)  and  Y  (translation),  of  the 
foil  are  resisted,  respectively,  by  linear  rotational  and  translation¬ 
al  spring  forces.  In  addition,  eddy  current  dampers  provide  forces 
for  increasing  the  damping  in  either  or  both  degrees  of  freedom. 
Signals  from  strain  gages  mounted  on  the  supporting  springs  give  a 
continuous  indication  of  displacements,  Y  and  6,  as  they  vary  with 
time;  see  Figure  1. 

I 

The  foil  itself  is  suspended  below  a  surface  plate  to  elimi¬ 
nate  wave  effects. 


FLUTTER  ANALYSIS 

Damping  ratio,  ‘  a  measure  of  the  rate  of  decay  of  oscilla¬ 
tions,  is  used  as  the  criterion  for  identifying  both  classical  and 
subcritical  flutter.  Analytic  expressions  for  damping  ratio  are  ob¬ 
tained  from  solutions  of  the  equations  of  motion  for  the  flutter  ap¬ 
paratus. 


EQUATIONS  OF  MOTION 

Two  separate  analyses  of  the  response  of  the  flutter  appa¬ 
ratus  --  the  Extended  Simplified  and  the  Modified  Theodorsen 
Analyses  --  were  proposed  by  McGoldrick  and  Jewell  in  Reference 
2. 


The  Modified  Theodorsen  Analysis  is  based  on  Theodorsen's 
solution  of  the  stability  problem  for  two-dimensional  flow  over  a 
plane  foil  of  infinite  aspect  ratio  with  no  structural  deunping.  ^ 
Modifications  of  the  equations,  originally  derived  by  Theodorsen, 
are  explained  fully  in  Reference  2.  Changes  were  made  in  expres¬ 
sions  for  the  lift  forces  cind  moments  amd  in  the  inclusion  of  a  steady 
hydrodynamic  moment  ascribed  to  circulation.  **  The  equations  of 

’^Damping  ratio  is  defined  as  the  ratio  of  "damping"  to  "critical 
damping"  (expressed  as  a  decimal)  for  an  oscillating  system. 

**Also  the  generally  complex  Theodorsen's  function,  C(k)  ,  was 
assumed  equal  to  1/2,  because  the  reduced  frequency  or  Strouhal 
number  (  <u  b/S  )  was  relatively  high  (of  the  order  of  1.  2).  The 
limited  calculations  carried  out  by  the  Hydromechanics  Laboratory 
showed  that  the  real  part  of  C(k)  was  very  nearly  equal  to  1/2 
and  that  the  phase  angle  of  C(k)  was  in  the  order  of  10  deg. 

These  findings  are  considered  to  be  in  agreement  vath  the  assumption. 
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motion  for  the  Modified  Theodors en  Analysis  are: 

1^*-,  (c  +  J-Ab2S)0  +  (k  -  ^AI.S2)fl- mhY  =  0 
2  2 

-mh0  -  AbS(9  -  iAS^fl  +  (M  +  m)  Y  +  (C  +  iAS)Y  +  KY  =  0 
2  2 

where  (')  and  (  '*  )  denote  first  and  second  total  derivatives  with 
respect  to  time.  * 

The  equations  of  motion  for  the  Extended  Simplified  Analy¬ 
sis  are: 

10  +  c0  +  (k  -  A1,S2)0  -  mhY  +  ALSY  =  0 

-mh0  -  AS^0  +  (M  +  in)Y  +  (C  +  AS)Y  +  KY  =  0 

For  both  the  Modified  Theodorsen  Analysis  and  the  Ex¬ 
tended  Simplified  Analysis,  the  equations  of  motion  are  mathemati¬ 
cally  of  the  form: 

Bij0  +  Bjj0  +  Bio0  +  BjjY  +  BjiY  =  0 

\\J  +  R„0  +  +  R,,Y  +  R,,Y  +  R,„Y  =  0 

For  this  reason,  the  following  methods  of  determining 
damping  ratios  may  be  applied  to  either  of  the  two  analyses  by  em¬ 
ploying  the  appropriate  coefficients. 

DIGITAL  COMPUTER  DETERMINATION  OF  DAMPING  RATIOS 

Assumed  solutions  of  the  form  0  =  and  Y  =  Y^e^^ 

(where  ^QfYg  and  \  are  generally  complex)  are  substituted  in  the 
equations  of  motion  for  both  the  Extended  Simplified  and  Modified 
Theodorsen  Analyses.  Two  algebraic  equations  of  the  form, 

(Bj2^^  +  ® lO^^O  ^  ^2 j A)Y Q  -  0 

(Rj2^*  +  Rj  jX  +  RjqI^q  +  (R22^^  ^20^^0  ~  ® 

result  from  this  substitution  for  either  analysis.  These  equations 
must  be  satisfied  so  that  the  assumed  solutions  may  satisfy  the  equa¬ 
tions  of  motion. 


’^AIl  symbols  in  the  coefficients  of  the  equations  of  motion  are  de¬ 
fined  in  the  notation,  p.  710. 
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A  nontrivial  solution  of  the  preceding  equations  exists  if, 
and  only  if,  the  determinant  of  the  coefficients  vanishes.  In  other 
words,  the  following  characteristic  equation  of  fourth  degree  (fre¬ 
quency  quartic)  in  the  generally  complex  frequency  A  =  ^  +  j<u  must 
be  satisfied: 


(Iij2^^22  ~  ^^11^^22  ^^12^21  ~  ^^11^22  “  ^12^21^^^ 

+  (Hjoll22  *^12*^0  “  *^10^22  ~  +  (BjqR2i  + 

+  R11H21  +  I^iill20  ~  *^0^21^^  *^10*^20  ~  ^10*^21  ~  ®  • 


Note  that  the  coefficients  of  this  equation  are  all  real. 

The  hydrofoil  of  the  flutter  apparatus  has  two  degrees  of 
freedom  and,  consequently,  two  modes  of  vibration  corresponding, 
in  general,  to  two  distinct  natural  frequencies.  If  one  of  the  two 
modes  is  considered  to  predominate  in  the  total  motional  response, 
the  damping  ratio  will  equal  the  logarithmic  decrement  of  the  re¬ 
sponse  curve  0  vs  t  and/or  Y  vs  t  divided  by  2n.  For  numerical 
values  of  damping  ratio  that  are  small  compared  to  unity,  damping 
ratio  may  also  be  approximated  by  -/x/w  .  ** 

Solutions  of  the  preceding  frequency  quartic  equation,  for 
a  given  set  of  coefficients,  may  occur  mathematically  as  four  real 

’^Logarithmic  decrement  S  may  be  considered  to  be  defined  by  the 
expression  S  =  2nc/C^  /yj\  -  (c/C^)*  •  When  c/C^  is  small  compared 

to  unity,  as  is  anticipated  in  the  present  study,  8  =  2nc/C^  .  Damp¬ 
ing  ratio,  defined  previously  as  c/C^,  is  then  approximately  equal 
to  the  logarithmic  decrement  8  divided  by  2n  . 

’^’’‘Consider  a  response  Y  =  Ype^‘  .  The  logarithmic  decrement  of 
the  Y  vs  t  curve  may  be  interpreted  as  the  natural  logarithm  of  two 
successive  maximum  amplitudes,  Yj  and  Y2,  which  occur  T  seconds 
apart.  "T"  is  therefore  the  period  of  oscillation;  i.  e.  T  =  2n/a)  , 
where  w  is  the  circular  frequency  of  the  oscillation.  Then 

5  I  ^0*^'  I  /i(ti-t2)  ,  -uT  TO/ 

R,(Y2)  Yoe^‘2 

Thus  damping  ratio,  which  is  equal  to  8/2n  ,  may  be  expressed  as 
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roots,  two  pairs  of  complex  conjugate  roots,  or  as  two  real  roots 
and  one  pair  of  complex  conjugate  roots. 

When  solutions  occur  as  two  pairs  of  complex  conjugates, 
two  of  the  roots  will  be  associated  with  one  mode  of  vibration  and 
the  remaining  two  will  be  associated  with  the  other  mode. 

For  oscillations  in  a  given  mode,  no  distinction  is  made  be¬ 
tween  damping  in  the  translational  degree  of  freedom  and  damping 
in  the  rotational  degree  of  freedom.  Thus  the  two  roots,  A  ,  asso¬ 
ciated  with  that  mode,  have  identical  values  of  and  w  ,  the  two 
roots  being  complex  conjugates. 

The  sign  of  the  damping  ratio,  -/i/w,  is  always  opposite  to 
the  sign  of  the  associated  value  of  damping  fi.  For,  a  negative 
value  of/i  ,  indicating  stability,  is  associated  with  a  positive  loga¬ 
rithmic  decrement  and,  hence,  a  positive  damping  ratio.  The  con¬ 
verse  statement  is  also  true. 

When  classical,  critical,  or  subcritical  flutter  occurs,  a 
marked  coupling  between  the  two  natural  modes  may  be  observed,  ^ 
depending  upon  the  initial  conditions.  For  the  two-degree- of- 
freedom  system,  this  is  mathematically  evident  from  a  considera¬ 
tion  of  the  four  distinct  eigenvalues  associated  with  the  equations  of 
motion  when  solutions  of  the  form  0  =  and  Y  =  Yoe^‘  are  as¬ 

sumed.  When  these  four  values  of  X  are  denoted  by  subscripts  1 
through  4,  the  total  response  may  be  expressed  as 


Y  —  V  V  0^2*^  .  Y  ^3^  1  V 

I  IgjC  +  I gjC  +  IgjB  +  1046  . 

If  any  of  the  Xj  (i  =  1,  2,  3,  4)  contains  a  positive  n  (un¬ 
stable  condition),  then,  regardless  of  the  initial  conditions  imposed 
upon  the  physical  system,  the  response,  0  vs  t  and  Y  vs  t  ,  must 
ultimately  be  unstable  for  the  coupled  mode  system.  That  is,  the 
unstable  terms  in  the  expression  for  0  and  Y  must  eventually  pre¬ 
dominate  in  the  total  response. 

The  frequency  quartic  equation  was  coded  for  solution  on 
an  IBM  704  electronic  data  processing  machine  by  the  TMB  Applied 
Mathematics  Laboratory.  A  sufficiently  wide  range  of  values,  for 
the  physical  parameters  that  mcdce  up  the  coefficients  of  the  fre¬ 
quency  quartic  equation,  was  employed  to  correspond  to 
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experimental  conditions  for  which  damping  ratio  vs.  speed  data  were 
available.  Solutions,  Aj  =  /i,  +  jw;  ,  of  the  frequency  quartic  equa¬ 
tion  were  obtained,  and  for  each  solution,  damping  ratio  and 

logarithmic  decrement  8  were  computed  as  part  of  the  program. 

The  results  of  these  digital  computations  are  presented  graphically 
in  Appendix  A.  Figure  2  shows  a  sample  format  of  the  output  data. 
The  corresponding  input  data,  together  with  the  coefficients  AQi  Ai, 
A2i  A3  of  the  frequency  quartic  equation,  have  also  been  included 
on  this  format.  Some  symbols  on  the  format  differ  from  the  notation 
used  throughout  this  report.  Table  1  illustrates  the  equivalence  of 
these  symbols. 

Note  that  the  results  of  the  computations  predict  the  sta¬ 
bility  or  instability  of  the  system  but  not  the  total  motional  response, 
because  Yg  and  are  not  found  if  the  initial  conditions  are  not  con¬ 
sidered  in  the  solution  of  the  equations  of  motion. 

ANALOG  DETERMINATION  OF  FLUTTER  RESPONSE 

The  equations  of  motion,  whether  based  on  the  Extended 
Simplified  Analysis  or  on  the  Modified  Theodorsen  Analysis,  may  be 
solved  by  an  electronic  differential  analyzer  to  obtain  a  set  of  com¬ 
pletely  predicted  response  curves,  6  =  9 it)  and  Y  =  Y(t).  Four  pos¬ 
sible  general  types  of  response  may  be  obtained;  see  Figure  3. 

These  include  stable  and  unstable  oscillations,  both  with  one  mode 
predominating  and  with  beating  between  the  two  modes  to  produce 
either  a  stable  or  an  unstable  condition.  "Temporarily  suppressed 
flutter,  "  a  special  case  of  the  previous  classification,  occurs  when 
the  damping  /x  is  positive  (unstable)  in  one  mode,  negative  (stable) 
in  the  other,  and  when  initial  conditions  are  such  as  to  med^e  the 
stable  mode  predominate  initially  in  the  total  response.  Eventually, 
the  unstable  mode  must  predominate  to  produce  classical  flutter. 
However,  the  number  of  cycles  which  must  pass  before  the  insta¬ 
bility  is  observed  will  depend  upon  the  physical  parameters  of  the 
apparatus.  Recognition  of  this  phenomenon  is  of  great  importance  in 
the  use  of  analog  methods  of  predicting  flutter  because  the  number 
of  cycles  over  which  the  response  curves  are  recorded  may  be  in¬ 
sufficient  to  allow  the  amplitudes  due  to  the  unstable  mode  to  reach 
observable  magnitudes.  The  difficulty  may  be  obviated  by  obtaining 
two  sets  of  response  curves  from  the  analog  for  each  set  of  para¬ 
meters,  the  first  with  initial  conditions  to  make  one  mode  predomi¬ 
nate  and  the  second  with  initial  conditions  to  make  the  other  mode 
predominate  in  the  total  response.  An  unstable  mode,  if  it  exists, 
is  then  easily  observed  after  a  few  cycles  from  either  of  the  response 
curves.  An  example  of  "temporarily  suppressed  flutter"  is  shown 
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TABLE  1 


Comparison  of  Notation  Used  in  Figure  2  with 
Notation  of  This  Report 


Notation  Defined 

on  Page  710  -  712 

Notation  of 
Figure  2 

C 

CC 

c 

C 

K 

KK 

k 

K 

M 

MM 

m 

M 

n/(i> 

MU/OMEGA 

MU 

0} 

OMEGA 

^0 

A-O 

^1 

A-1 

Az 

A- 2 

As 

A- 3 

mh 

MH 

L 

L 

S 

S 
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Time 

SUble«One  Mode  Predominating 
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Unetable-one  Mode  Predominating 


Time 

Stabte-Beating  Between  Two  Modes 


■vw>/\/W/vi''y\/v 

j 

Tim* 

UniUbl*-B**tlnf  iMtwMn  Two  Mod** 


Temporarily  Suppraaaed  Fluttar 


Figure  3  -  Analog  Response 
Showing  Stable  and  Unstable 
Oscillations 
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in  Figure  3. 

In  the  present  study,  the  analog  approach  was  used  to  ob¬ 
tain  qualitative  information  about  the  general  character  of  the  re¬ 
sponse  wave  form.  However,  numerical  values  of  damping  ratio 
may  be  computed  from  response  curves,  in  which  one  mode  pre¬ 
dominates,  by  evaluating  b/2p  ,  where  B  is  the  logarithmic  decre¬ 
ment  of  these  curves. 

EXPERIMENTAL  WORK 

The  earliest  data  obtained  from  experiments  with  the  TMB 
Control  Surface  Flutter  Apparatus  were  presented  in  Reference  2. 
The  damping  ratios  and  associated  frequencies  therein  were  com¬ 
puted  from  records  of  vibrations  induced  in  the  apparatus  by  striking 
it  laterally  while  underway.  In  these  tests,  the  vertical  axis  of  rota¬ 
tion  was  located  at  the  forward  quarter-chord  position. 

A  more  detailed  explanation  of  the  experimental  procedure 
can  be  found  in  Reference  2. 

Since  the  publication  of  Reference  2,  the  Hydromechanics 
Laboratory  has  conducted  further  tests  on  the  flutter  apparatus.  The 
results  of  those  tests  have  been  published  in  Reference  6  and  are  in¬ 
cluded  in  this  study  for  comparison  with  the  results  predicted  by  the 
Extended  Simplified  and  Modified  Theodorsen  Analyses.  Table  2 
summarizes  the  various  configurations  for  which  the  comparison 
was  made.  Each  configuration  corresponds  to  the  different  values 
of  the  parameters  as  shown.  The  most  significant  changes  in  these 
parameters  were  an  increase  in  the  mass  of  the  assembly  that  ro¬ 
tates  with  the  hydrofoil  and  the  relocation  of  the  center  of  mass  fur¬ 
ther  aft  of  the  vertical  axis  of  rotation. 

Classical  flutter  was  observed  in  tests  of  Configurations  B 
and  D.  Experimentally  observed  damping  ratios  for  these  configura¬ 
tions,  as  well  as  for  the  configuration  (TMB  1222)  used  in  the  study 
reported  in  Reference  2,  are  compared  with  auia],ytical  predictions 
in  the  following  section. 

AhfALYTICAL  AND  EXPERIMENTAL  RESULTS 

Computations  of  damping  ratio  and  frequency  were  made 
for  Configurations  B,  D,  and  TMB  1222,  based  on  the  Extended 
Simplified  Analysis  and  the  Modified  Theodorsen  Analysis.  Numer¬ 
ical  values  of  the  parameters  mh,  L,  S,  and  a^,  for  which  the 
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computations  were  made,  are  given  in  Appendix  B.  Figures  19 
through  86,  presented  in  Appendix  A,  are  plots  of  computed  damping 
ratio  and  associated  frequency  <u  vs.  speed  for  both  natural  modes. 

Empirical  data  published  in  Reference  2  and  more  recent 
data  furnished  by  the  TMB  Hydromechanics  Laboratory^  are  super¬ 
posed  upon  these  plots  to  permit  the  desired  comparison  between 
analytical  predictions  and  experimental  observations  of  stability. 

The  most  significant  information  contained  in  these  graphs, 
for  purposes  of  comparison  of  theory  with  experiment,  has  been  re¬ 
produced  in  Figures  4  through  16.  The  section.  Discussion  of 
Results,  gives  the  reasons  for  selecting  these  graphs. 

Figures  17  and  18  show  the  experimentally  observed  de¬ 
pendence  of  critical  flutter  speed  on  mass  unbalance,  mh  for  Con¬ 
figurations  B  and  O,  respectively.  Superposed  on  those  figures  are 
analytical  curves  of  critical  flutter  speed  vs.  mh,  based  on  the 
Modified  Theodorsen  Analysis.  Note  that  these  curves  are  based  on 
different  values  of  L  even  though  ao  =1.  88^  for  all  tests  conducted 
for  Configurations  B  and  D.  The  use  of  values  of  L  not  uniquely  de¬ 
pendent  on  Aq  violates  the  assumption^  which  underlies  the  Modi¬ 
fied  Theodorsen  and  Extended  Simplified  Analyses;  see  Discussion 
of  Results  and  Recommendations. 

DISCUSSION  OF  RESULTS 

The  data  in  Figures  4  through  16  have  been  selected  from 
Figures  19  through  86  in  Appendix  A  for  convenience  in  interpreting 
the  comparison  of  analytical  and  experimental  damping  ratios.  The 
following  points  form  the  criteria  upon  which  the  data  of  Figures  4 
through  16  were  selected: 

1.  Not  all  values  of  L,  for  which  families  of  damping  ratio 
and  frequency  vs.  speed  curves  were  drawn  in  Appendix  A,  corre¬ 
spond  to  experimental  conditions.  For  Configuration  TMB  1222, 
the  observed^  value  of  L  was  2.  8  in. ,  corresponding  to  «o  =  0^. 

No  experimental  measurement  of  L  was  obtained  for  Configurations 
B  and  D.  Therefore,  based  on  computations,  a  family  of  curves  for 
the  parameter  L  was  plotted  in  am  attempt  to  find  a  correlation  be¬ 
tween  analytical  curves  for  a  unique  value  of  L  and  experimental 
data.  This  approach  for  Configurations  B  and  D  is  discussed  in 
Point  3. 


2.  For  Configuration  TMB  1222,  where  no  classical 
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Figure  4  -  Comparison  of  Computed  and  Experimental  Damping 
Ratios  (Configuration  TMB  1Z22,  mh  =  0.0  lb-sec^) 


728 


Damping  Rjtto  Dftmpint  Ratio 


Leibowitz  and  Belz 


Confifuratton  TMB*1222  mh*0.3  Ao'O  Lower  Mode 


Figure  5  -  Comparison  of  Computed  and  Experimental  Damping 
Ratios  (Configuration  TMB12Z2,  mh  =  0.3  lb-sec  ) 
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Figure  6  -  Comparison  of  Computed  and  Experimental  Damping 
Ratios  (Configuration  TMBlZ22,mh  =  2.0  lb-sec^) 


3 


Leibowitz  and  Belz 


ConftKuratlon  B  mh*3.53  Hitfhtr  Mod* 


Figure  7  -  Comparison  of  Cozppbted  and  Experimental  Damping 
Ratios  (Configuration  B,  mh  =  3.53  lb-sec^) 
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Figure  8  -  Comparison  of  Computed  and  Experimental  Damping 
Ratios  (Configuration  B,  mh  =  3.68  lb-sec^) 
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Figure  9  -  Comparison  of  Computed  and  Experimental  Damping 
Ratios  (Configuration  B,  mh  =  3,82  lb- sec  ) 
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Figure  10  -  Comparison  of  Computed  and  Experimental  Damping 
Ratios  (Configuration  B,  mh  =  3.94  lb-sec^) 
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Figure  11  -  Comparison  of  Computed  and  Experimental  Damping 
Ratios  (Configuration  B,  mh  =  4.06  lb-sec^) 
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Figure  12  -  Comparison  of  Computed  and  Experimental  Damping 
Ratios  (Configuration  B,  mh  =  4.21  lb-sec^) 
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Figure  13  -  Comparison  of  Computed  and  Experimental  Damping 
Ratios  (Configuration  B,  mh  =  4.32  lb-sec^) 
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Figure  14  -  Comparison  of  Computed  and  Experimental  Damping 
Ratios  (Configuration  D,  mh  =  4.05  lb-sec^) 
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CondfEurttlon  D  mh  •  4.21  Htflitr  Modt 


Figure  15  -  Comparison  of  Computed  and  Experimental  Damping 
Ratios  (Configuration  D,  mh  =  4.21  lb-sec^) 
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Figure  16  -  Comparison  of  Computed  and  Experimental  Damping 
Ratios  (Configuration  D,  mh  =  4.32  lb-sec^) 
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Figure  17  -  Comparison  of  Computed  and 
Experimental  Critical  Flutter  Speeds  for  an 
Assumed  Value  of  L  (Configuration  B) 
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Figure  18  -  Comparison  of  Computed  and 
Experimental  Critical  Flutter  Speeds  for  an 
Assumed  Value  of  L  (Configuration  D) 


742 


Leibowitz  and  Belz 


flutter  was  observe,  the  phase  angle  between  Y  and  0  indicates  that 
the  mode  of  lower  frequency  predominated  in  the  total  response  for 
mh  =  0.  0  and  0.  3  lb-sec^;  however,  for  mh  =  2.  0  lb-sec^  the  mode 
of  higher  frequency  predominated;  see  Appendix  C.  Therefore,  com¬ 
parisons  of  theory  and  experiment  should  be  made  for  the  lower  mode 
when  mh  =  0.  0  or  0.  3  lb-sec^  and  fo#  the  higher  mode  when  mh  = 

2.  0  lb-sec^’ 

3.  In  Configurations  B  and  D,  for  all  values  of  mh  under 
consideration,  the  higher  mode  consistently  gave  the  best  compari¬ 
son  with  test  data  and  is  therefore  postulated  to  have  been  the  mode 
predominantly  excited  in  the  total  response. 

4.  In  Figures  4  through  16,  damping  ratio  vs.  speed  is 
plotted;  Figures  19  through  86  in  Appendix  A  show  frequency  vs. 
speed  in  addition.  These  frequency  curves ‘were  omitted  in  Figures 
4  through  16  not  because  they  are  irrelevant  but  because  the  major 
concern  in  plotting  the  latter  figures  was  to  effect  a  comparison  of 
stability  predictions  (damping  ratios)  and  observations.  A  compari¬ 
son  of  analytical  frequency  predictions  with  observed  frequencies 
forms  an  important  part  of  the  general  comparison  of  theory  with 
experiment,  but  does  not  directly  offer  a  comparison  of  stability 
predictions  and  observations. 

Figures  4  through  6  furnish  convincing  evidence  that  re¬ 
sults,  based  on  the  Modified  Theodorsen  Analysis,  compare  well 
with  test  data  for  the  lower  range  of  values  of  mass  unbalance  (mh) 
under  consideration.  *  For  the  two  lowest  values  of  mh  (Figures  4 
and  5),  the  departure  of  analytic  damping  ratio  curves  from  experi¬ 
mental  data  is  generally  less  than  25  percent  of  the  experimental 
values.  At  mh  =  2.  0  lb-sec^,  the  scatter  of  test  data  is  much  great¬ 
er  than  at  the  lower  mh  values;  yet  the  analytical  curve,  based  on 
the  Modified  Theodorsen  Analysis,  predicts  damping  ratios  near  the 
average  of  scattered  empirical  data. 

In  the  consideration  of  Figures  31  through  86,  from  which 
Figures  7  through  16  have  been  abstracted,  no  unique  value  of  L  pro¬ 
vides  analytical  curves  that  correlate  with  test  data  for  Configura¬ 
tions  B  and  D.  This  lack  of  correlation  is  associated  with  values  of 
mass  unbalance  (mh)  about  one  and  one-half  to  two  times  as  large 
as  the  mh  values  of  Configuration  TMB  1222,  for  which  good  corre¬ 
lation  was  obtained  between  test  data  and  the  Modified  Theodorsen 

*That  is,  for  mh  =  0.  0,  0.  3,  and  2.  0  lb-sec^  (Configuration 
TMB  1222). 
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Analysis.  Because  experimental  estimates  of  L  were,  unfortunately, 
not  obtained  for  Configurations  B  and  D,  the  correlation  considered 
here  is  that  between  the  region  in  which  test  data  fall  and  a  "pre¬ 
dicted  region,  "  considered  to  be  bounded  by  two  curves  of  the  family 
of  curves  for  the  parameter  L.  An  arbitrary  choice  of  -1  and  -4 
as  values  of  L  provides  the  bounding  curves  in  which  test  data  gen¬ 
erally  fall;  see  curves  shown  in  Figures  7  through  16.  These  bounds 
are  associated  with  a  rather  narrow  range  of  variation  in  damping 
ratio;  the  differences  in  ordinates  of  the  damping  ratio  vs.  speed 
curves  based  on  the  Modified  Theodorsen  Analysis  for  L  =  -1  and 
L  =  -4  rarely  exceed  0.  02  at  a  given  speed  below  15  knots.  For  ex¬ 
ample,  in  Figure  13,  the  Modified  Theodorsen  Analysis  predicts  a 
damping  ratio  of  0.  004  at  a  speed  of  12.  9  knots  for  L  =  -1  and  a 
damping  ratio  of  /-O.  014  for  L  =  -4  at  that  speed.  The  difference  in 
the  ordinates  for  this  example  is  thus  0.  018. 

Reference  2  reported  that  at  a  preset  attack  angle  Og  ** 
of  zero  degrees,  both  steady  lift  forces  and  moments  were  recorded. 
One  important  feature  of  both  the  Modified  Theodorsen  Analysis  and 
the  Extended  Simplified  Analysis  is  that  they  include  a  term  for 
steady  hydrodynamic  moment  due  to  circulation  to  account  for  this 
experimental  observation  on  the  TMB  Flutter  Apparatus  with  the  axis 
of  the  rudder  stock  at  the  forward  quarter-chord  position.  The  clas¬ 
sical  Theodorsen  analysis^  predicts  zero  steady  hydrodynamic  mo¬ 
ment,  i.  e.  ,  hydrodynamic  balance  when  the  stock  is  at  this  position, 
thus  disagreeing  with  experimental  observation.  The  presence  or 
absence  of  hydrodynamic  balance  makes  a  marked  difference^  in  the 
overall  damping  characteristic  obtained  from  the  flutter  apparatus. 

It  should  be  recognized  that  the  hydrofoil  of  the  flutter  apparatus  has 
an  NACA  0015  cross  section,  i.  e.  ,  its  thickness  is  15  percent  of 
the  chord  length,  which  may  well  be  a  significant  departure  from 
the  plane  foil  assumed  in  the  classical  theory. 

Nevertheless,  it  has  not  been  conclusively  established 
that  the  preset  attack  angle  was  actually  zero,  because  mecheui- 
ical  imperfections  of  the  flutter  apparatus  auid  its  mounting  intro¬ 
duce  uncertainties  as  high  as  0.  5  deg  in  nominal  values  of  .  The 

’’‘The  parameter  L,  the  moment  arm  of  the  steady  lift  force  about 
the  axis  of  rotation,  was  defined  as  the  ratio  of  observed  steady 
moment  to  observed  steady  lift  force. 

**  The  preset  attack  angle  defines  the  equilibrium  orientation  of 
the  hydrofoil;  see  Figure  1. 
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values  of  og  considered  in  this  report  are  small  {  og  =  0.0  deg. 
and  1.  88  deg).  Therefore,  uncertainties  in  must  be  considered 
as  a  possible  explanation  of  the  apparently  observed  nonzero  lift  and 
moment  at  Og  =  0°. 

There  is  no  conclusive  evidence,  either  from  experiments 
or  theoretical  considerations,  to  cause  one  or  the  other  of  the  pre¬ 
ceding  explanations  to  account  for  the  observation  of  lift  forces  and 
moments  at  Oq  =  0  deg,  as  reported  in  Reference  2.  For  further 
discussion  of  the  parameter  L  associated  with  these  moments,  see 
Appendix  D. 

When  applied  to  the  TMB  Flutter  Apparatus,  the  classical 
Theodorsen  equations:^'  ^ 

1.  Yield  a  steady  lift  force  about  three  times  as  large  as 
that  found  experimentally:  i.  e.  ,  the  computed  force  and  moment  do 
not  converge  to  the  value  of  the  experimental  steady  lift  force  and 
moment  at  zero  frequency.  Hence,  they  cannot  predict  divergence, 
which  is  usually  included  in  a  flutter  analysis. 

2.  Show  no  dependence  on  the  mean  angle  of  attack,  where¬ 
as  the  steady  moment  was  found  to  vary  with  the  mean  angle  of 
attack. 


3.  Show  good  agreement  with  observed  critical  speeds 
over  a  limited  range  of  mh;  see  Reference  7.  However,  this  refer¬ 
ence  gives  no  information  on  agreement  of  predicted  frequencies  auid 
phase  relations.  A  sound  method  should  check  in  all  respects. 

It  should  be  noted  that  the  majority  of  the  conditions  ex¬ 
plored  experimentally  and  reported  in  Reference  7  show  the  "mild 
flutter"  condition  in  which  the  inclusion  or  omission  of  structural 
damping  caui  make  a  relatively  large  difference  in  the  predicted  criti¬ 
cal  flutter  speed. 

Analytical  curves  of  damping  ratio  vs.  speed  based  on 
the  Extended  Simplified  Analysis  do  not  fall  along,  or  generally  near, 
the  experimental  data  plotted  except  for  the  two  lowest  values  of  mass 
unbalance  (mh)  employed;  i.  e.  ,  mh  =  0.  0  and  0.  3  lb-sec^;  see  Fig¬ 
ures  4  through  6. 

For  analog  solutions  of  the  equations  of  motion,  initial 
conditions  imposed  on  the  foil  must  be  taken  into  consideration.  In 
general,  when  arbitrary  initial  conditions  are  assumed,  the  solutions 
Y  (t)  and  0(i)  of  the  differential  equations  of  motion  for  the 
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flutter  apparatus  represent  the  superposition  of  two  modes  of  vibra¬ 
tion,  The  solution  can,  however,  be  restricted  to  represent  a 
single  mode  of  vibration;  one  of  the  two  modes  can  be  suppressed,  if 
the  initial  conditions  are  properly  chosen.  This  will  eliminate  the 
possibility  of  beating  between  modes,  which  is  the  cause  of  "tempo¬ 
rarily  suppressed  flutter,  "  discussed  under  Analog  Determination  of 
Flutter  Response.  '''  Criteria  for  the  proper  selection  of  initial  con¬ 
ditions  to  cause  one  mode  to  predominate  in  the  total  response  are 
established  in  Appendix  E.  For  a  discussion  of  initial  conditions  as 
they  affect  the  vibration  of  a  ship's  rudder,  see  Appendix  F. 

In  Figures  17  and  18  the  regions  in  which  experimental 
and  analytical  data  fall  are  quite  close  and  sometimes  overlap. 

This  suggests  a  strong  dependence  of  critical  flutter  speed  upon  the 
parameters  mh  and  L.  Note  that  a  good  correlation  of  observed  and 
analytical  critical  speeds  as  functions  of  mh  are  possible  only  by 
violating  the  assumption,  made  in  Reference  2,  that  L  is  a  function 
of  Oq  only.  For,  in  Figure  17  a  value  of  L  =  -3.  0  was  assumed  to 
correspond  to  Oq  =  1.  88  deg,  whereas  in  Figure  18,  a  value  of 
Li  =  -1.  5  was  assumed  to  correspond  to  a,,  =  1.  88  deg.  Therefore, 
the  "correlation"  between  theory  and  experiment  shown  in  Figures 
17  and  18  has  no  meaning  other  than  to  indicate  that  L  is  apparently 
a  function  of  mh  as  well  as  of  .  (see  Recommendation  2). 

CONCLUSIONS 

The  preceding  comparison  of  analytical  stability  predic¬ 
tions  with  available  test  data,  although  quite  limited  in  extent,  never¬ 
theless  leads  to  the  following  conclusions: 

1,  The  Modified  Theodorsen  Analysis,  as  compared  with 
the  Extended  Simplified  Analysis,  appears  to  be  the  most  suitable 
analysis  for  yielding  good  predictions  of  damping  ratio  and  frequency 
for  a  given  speed  amd  for  predicting  critical  flutter  speeds. 

2.  The  Extended  Simplified  Analysis  correlates  well  with 
test  data  for  only  the  two  lowest  values  of  mh  investigated;  i.  e.  , 
mh  =  0.  0  and  0.  3  lb-sec^.  This  amalysis  is  therefore  considered  to 
be  quite  limited  in  its  applicability  to  the  problem  of  predicting 
flutter. 


3,  The  accuracy  of  the  Modified  Theodorsen  Analysis 

’“‘This  does  not  imply  that  "temporarily  suppressed  flutter"  will 
necessarily  occur  when  two  modes  beat. 
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apparently  improves  as  mass  unbalance  mb  is  decreased.  The  de¬ 
pendence  of  L  on  parameters  other  than  (such  as  mh  or  speed) 
might  be  profitably  investigated  to  extend  the  range  of  parameters 
over  which  the  Modified  Theodorsen  Analysis  correlates  well  with 
test  data  (see  Recommendation  2). 

4.  The  ship  designer  requires  the  simplest  analysis  that 
will  yield  predictions  consistent  with  experimental  observations.  The 
results  obtained  for  theoretical  and  experimental  data  indicate  that 
the  Modified  Theodorsen  Analysis  might  yield  more  realistic  pre¬ 
dictions  for  ship  control- surface  systems  than  the  more  elaborate 
equations  used  in  the  classical  flutter  theory. 

RECOMMENDATIONS 

To  further  the  development  of  practical  design  criteria 
for  avoiding  control- surface  flutter,  the  following  recommendations 
are  made: 


1.  A  series  of  hydrofoil  towing  tests  should  be  conducted 
to  determine 

(1)  the  point  of  application  of  the  hydrodynamic  lift 
force  on  the  foil  and 

(2)  whether  steady  lift  and  moment  do  occur  at 

Og  =  0  deg.  These  tests  may  not  be  possible  on  the  existing  TMB 
Control  Surface  Flutter  Apparatus  because  of  the  uncertainties  in 
measurements  of  Og  that  are  inherent  in  the  present  device.  There¬ 
fore,  in  the  design  of  future  hydrofoil  flutter  apparatus,  considera¬ 
tion  should  be  given  to  methods  of  measuring  the  chordwise  location 
of  center  of  pressure  for  variations  in  og  ;  see  recommendation  2. 

It  is  recommended  at  this  stage  that  the  value  of  L,  the  distance  be¬ 
tween  the  axis  of  rotation  and  the  center  of  pressure,  be  determined  by 
taking  the  ratio  of  the  observed  moment  and  lift  force  as  indicated 
in  Reference  2. 

2.  If  the  location  of  the  lift  force  as  determined  in  accord¬ 
ance  with  Recommendation  1,  is  found  to  depart  significantly  from 
the  forward  quarter-chord  position,  a  series  of  tests  should  be  car¬ 
ried  out  to  determine  the  dependence  of  L  on  other  parameters,  as 
well  as  on  Og  .  The  result  of  this  recommended  study  will  be  a 
graphical,  or  possibly  an  analytical,  relation  between  L  and  these 
parameters.  In  addition,  the  assumption  that  the  oscillatory  lift 

and  moment  follow  the  steady  relations  should  be  investigated.  The 
more  realistic  expression  for  L  so  obtained  might  then  be  used  in 
conjunction  with  the  Modified  Theodorsen  Analysis  to  extend  the 
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range  of  physical  parameters  over  which  that  analysis  agrees  well 
with  experimental  findings. 

3.  A  study  of  the  applicability  of  the  classical  Theodorsen 
approach  to  flutter  prediction  for  model  and  full-scale  rudders 
should  be  undertaken  along  with  a  continuing  study  of  both  the 
Extended  Simplified  and  Modified  Theodorsen  Analyses.  Investiga¬ 
tion  of  the  latter  two  analyses  is  particularly  stressed  because  of 
the  relative  ease  with  which  they  may  be  applied  and  because  they 
are  formulated  in  dimensional  terms  common  to  present  vibrations 
analysis  practice.  It  would  be  of  value  to  apply  the  classical  Theo- 
dorson  equations  to  the  experimental  data  published  in  this  report. 

A  model  devised  for  flutter  investigation  should  be  designed 
so  that  its  experimental  operating  conditions  and  (certain)  parameter 
values  are  similar  to  those  for  the  full-scale  rudder;  see  Appendix 
G. 


In  evaluating  the  applicability  of  the  classical  Theodorsen 
or  other  analyses  for  flutter  prediction,  it  is  desirable  to  experi¬ 
mentally  determine  the  reaction  on  vibrating  rudders  and/or  rudder 
models.  (In  Appendix  H  an  analysis  is  given  for  the  reduction  of 
the  classical  Theodorsen  results  to  concrete  form.  The  dependence 
of  the  derived  formulas  on  ship  parameters  is  shown,  and  a  discus¬ 
sion  is  given  of  the  applicability  of  Theodorsen's  results  to  rudders 
and  the  need  for  measurements  of  the  reactions  on  the  rudder). 

4.  The  discrepancy  between  stability  predictions  based  on 
the  Routh  discriminant’^  of  the  frequency  quartic  equation,  and  on 
solutions  of  that  same  equation  as  reported  in  the  present  study, 
should  be  resolved. 

5.  The  dependence  of  critical  flutter  speed  on  mass  un¬ 
balance  should  be  investigated. 

6.  It  is  suggested  that  dimensional  notation  be  used  in  any 
hydroelastic  analysis  made  for  consistency  with  existing  codes  at 
the  Applied  Mathematics  Laboratory  on  hull  vibration  analysis. 

While  dimensionless  notation  has  found  wide  acceptance  by  naval 

’^Stability  predictions  for  Configuration  TMB  1222,  based  on  the 
Routh  discriminant,  can  be  found  in  Reference  2,  where  a  discrep¬ 
ancy  is  also  shown  to  exist  between  the  Routh  criteria  for  stability 
and  experimental  results. 
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architects  and  aeronautical  engineers,  it  is  not  widely  used  in  the 
analysis  of  mechanical  vibration  at  the  present  time. 
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APPENDIX  A 
COMPARISON  OF 

PREDICTED  AND  OBSERVED  DAMPING  RATIOS 
AND  ASSOCIATED  FREQUENCIES 
FOR  THE 

TMB  CONTROL  SURFACE  FLUTTER  APPARATUS 


Analytical  curves  of  damping  ratio  and  frequency  vs.  speed, 
obtained  from  digital  solutions  of  the  frequency  quartic  equation 
based  on  both  the  Extended  Simplified  and  Modified  Theodorsen 
Analyses,  are  ^presented  in  Figures  19  through  86  for  combinations 
of  the  parameters  given  in  Table  3. 

Superposed  on  these  plots  are  experimental  values  of  damp¬ 
ing  ratio  and  frequency  vs.  speed.  The  comparison  of  theory  with 
experiment  should  be  made  in  accordance  with  the  points  given  in 
the  Discussion  of  Results. 

Where  experimental  data  were  not  available  (Figures  59 
through  74),  analytical  predictions  are  given  to  present  a  qualita¬ 
tive  picture  of  trends,  in  the  computed  curves,  produced  by  changes 
in  the  parameters  mh,  L,  and  S.  (Note  that  in  the  present  analy¬ 
ses,  L  is  a  function  of  Og  only.  ) 
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TABLE  3 

Numerical  Values  of  Parameters  mh,  L,  S,  and 

for  which 

Analytical  Frequency  and  Stability  Predictions 
Have  Been  Made 


Parameters 

Dimensions 

Numerical  Values 

Configuration  TMB 

1222 

mh 

lb-sec^ 

0.  0,  0,  3,  2.  0 

L 

inches 

0.  0,  2.  8,  6.  41,  8.  05, 

8.  90,  10.  0 

S 

knots 

0.  0  to  20.  0  in  1-knot 
increments 

«o 

degrees 

0.  0 

Configurations  B  e 

,nd  D 

mh 

lb-sec^ 

3.  53,  3.  68,  3.  82,  3.  94, 

4.  05,  4.  21,  4.  32 

L 

inches 

-4,  -3,  -2,  -1,  0,  2.8, 

6.  41,  8.  05,  8.  90,  10.  0 

S 

knots* 

0.  0  to  20.  0  in  1-knot 
increments 

«o 

degrees 

1.  88 

In  performing  calculations,  S  must  be  in  units  consistent  with 
those  of  the  other  parameters;  i.  e.  ,  in  inches  per  second.  For 
convenience  of  the  reader,  speeds  are  indicated  in  knots  in 
Figures  4  through  86. 

For  each  configuration,  all  possible  permutations  of  the  above 
parameters  were  considered  analytically. 
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Figure  21 
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Figure  23 
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Figure  24 
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Figure  79 
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Figure  80 
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Figure  84 
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Figure  85 
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Figure  86 
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APPENDIX  B 

NUMERICAL  VALUES  OF 
THE  PARAMETERS  mh,  L,  S,  and 

The  numerical  values  of  mh,  L,  S,  and  «o  >  given  in 
Table  3,  are  those  values  for  which  analytical  stability  predictions 
were  made,  based  on  the  Extended  Simplified  and  Modified  Theodor- 
sen  Analyses. 

The  mass  unbalance  (mh)  values  of  Configurations  B,  D,  and 
TMB  1222  were  chosen  to  correspond  to  experimental  values  of  mh. 

For  mh  values  used  in  conjunction  with  Configuration  TMB 
1222,  see  Reference  2.  For  Configurations  B  and  D,  mh  values 
were  computed  from  numerical  values  of  ,  m'  ,  and  ab 
supplied  by  the  Hydromechanics  Laboratory.  The  relation  between 
these  quantities  is 


mb  =  -  >»*  (abK 

Where  ab  is  the  location  of  the  axis  of  rotation  aft  of  the  midchord 
position,  m'  is  the  experimentally  observed  added  mass  in  transla¬ 
tion,  and  is  the  moment  about  the  forward  quarter-chord  point 
of  the  rotational  structural  mass.  For  both  Configurations  B  and  D, 
ab  =  -4.  5  in.  and  m'  =  0.  284  Ib-sec^/in.  Values  of  S„  employed 
experimentally  were  2.  25,  2.  40,  2.  54,  2.  66,  2.77,  2.  93,  and  3.  04 
lb-sec^. 


The  values  of  L  listed  under  Configuration  TMB  1222  are 
all  those  values  for  which  stability  computations  based  on  the  Routh 
discriminant  were  made  in  Reference  2.  Experimental  measure¬ 
ments  reported  in  that  reference  showed  that  L  =  2.  8  in.  was  the 
experimental  value  to  be  considered  for  ao  =  0  deg. 

No  experimental  measurements  of  L  were  obtained  for  Con¬ 
figurations  B  and  D.  A  range  of  L  corresponding  to  that  employed 
for  Configuration  TMB  1222  was  used  in  addition  to- the  negative 
values,  L  =  -1,  -2,  -3,  -4,  which  were  arbitrarily  selected  to  ex¬ 
amine  the  trends  of  computed  damping  ratio  and  frequency  vs.  speed 
curves. 


The  range  of  speeds  S  corresponds  to  the  design  range  of 
speeds  of  the  TMB  Control  Surface  Flutter  Apparatus  (0  to  20  knots). 


820 


Leibowitz  and  Belz 


The  no’**  listed  in  Table  3  are  equal  to  nominal  set¬ 
tings  of  «o  which  tests  were  conducted  on  the  flutter  apparatus. 


APPENDIX  C 


RELATIONSHIP  OF  THE  MODAL  FREQUENCIES 
TO  THE  PHASE  ANGLE  BETWEEN  Y  and  0 

In  the  ideal  two-degree-of-freedom  system  shown  in  Fig¬ 
ure  87A,  the  fundamental  or  lower  frequency  is  observed  when  the 
motions  of  the  two  masses  are  in  phase.  ®  When  the  motions  are  in 
phase,  the  algebraic  signs  of  and  X2  are  the  same  at  any  instant 
of  time. 


When  analyzing  the  motions  of  the  flutter  apparatus,  the 
fundamental  frequency  is  also  observed  when  Y  and  0  have  the 
same  signs.  However,  this  condition,  for  the  flutter  apparatus, 
indicates  a  phase  angle  <!>  of  180  deg. 

This  situation  is  explained  by  considering  the  sign  conven¬ 
tions  employed  in  the  two  systems  shown  in  Figure  87.  Table  4 
shows  the  relation  between  frequency  and  phase  angle  for  both  sys¬ 
tems. 


Experimental  values  of  phase  angle  <i>  and  frequency  are 
given  in  Table  4  of  Reference  2  for  Configuration  TMB  1222. 
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AT  REST 


IN  PHASE  (<A  =  0°) 


OUT-OF- PHASE  (6  =  180°) 


Figure  87a  -  Two-Mass,  Two-Degree-of-Freedom  System 


AT  REST  IN  PHASE  (0  =  0°)  OUT  OF  PHASE  (<f>  =  180°) 

Figure  87b  -  Hydrofoil  of  TMB  Control  Surface  Flutter  Apparatus 


Figure  87  -  Sign  Convention  for  Two-Degree-of-Freedom 
Mass  Spring  System  and  Hydrofoil  System 
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TABLE  4 

Relation  between  Frequency  and  Phase  Angle  (0°) 

for 

Connmon  Two- Mass  System  and 
Hydrofoil  of  Flutter  Apparatus 


Two- Mass  System 
(Figure  87a) 

Flutter  Apparatus  Hydrofoil 
(Figure  87b) 

6  Deg 

Signs  of 

Xj  and  Xj 

Frequency 

(f)  Deg. 

(angle  between 
Y  and  $  ) 

Signs  of 
Y  and  6 

Fre¬ 

quency 

m 

Same 

Fundamental 

(Lower) 

=  0° 

Same 

Higher 

D 

Opposite 

Higher 

=  180° 

Oppos. 

Funda¬ 

mental 
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APPENDIX  D 


RELATIONSHIP  BETWEEN  THE  PARAMETER  L 
AND  THE  ANALYTICAL  EXPRESSIONS 
FOR  HYDRODYNAMIC  LIFT  AND  MOMENT 


During  constant  speed  runs  of  the  flutter  apparatus  with  the 
hydrofoil  at  a  nominal  zero  preset  angle  of  attack,  steady  lift  forces 
and  moments  were  observed  by  the  authors  of  Reference  2.  ’  To 
account  for  this  observation  it  was  assumed  that  the  center  of  lift 
did  not  coincide  with  the  axis  of  rotation  at  the  forward  quarter-chord 
position.  2  The  observation  of  steady  forces  and  moments  cannot  be 
accounted  for  by  the  aerodynamic  strip-foil  theory  used  to  design 
the  flutter  apparatus.  Also,  the  assumption  that  the  center  of  lift 
did  not  coincide  with  the  forward  quarter-chord  position  was  made 
contrary  to  strip  theory. 

For  both  the  Extended  Simplified  and  Modified  Theodors en 
Analyses,  it  was  assumed  that  the  component  of  oscillatory  moment 
Mg  ,  due  to  the  oscillatory  lift  force,  follows  the  steady  moment 
relation;  i.  e, ,  Mg  ~  . 

For  the  Extended  Simplified  Analysis,  the  oscillatory  lift 
force  is  taken  to  be 


which  is  an  expression  for  the  instantaneous  oscillatory  lift  force. 
The  total  oscillatory  moment  was  considered  to  be  caused  by  F,^ 
acting  at  some  point  called  the  center  of  pressure,  not  on  the  axis 
of  rotation.  In  other  words,  for  the  Extended  Simplified  Analysis, 

Mg  =  LFl 

or 

Mg  =  ALS^^0-I^. 

According  to  the  classical  Theodorsen  analysis,  as  applied 
to  the  flutter  apparatus,  ^  the  oscillatory  lift  forces  and  moments 
should  be  described  by  the  expressions 

A  great  deal  of  the  material  contained  in  this  appendix  is  a  restate¬ 
ment  of  portions  of  Reference  2.  This  material  is  included  here  for 
the  reader's  convenience. 
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and 


K  i,\s2rt-i  \SV  +  \bSfl 
I.  2  2 

Mfl  =  -  T  '‘‘'S'’ 


The  Modified  Theodorsen  Analysis  retains  the  classical 
expression  for  F^.  But  in  describing  ,  both  the  classical  ex¬ 

pression  (-  -  \l)‘SiO)  and  a  component  to  account  for  the  moment 
due  to  are  considered.  Thus 

\1^  =  -  1  +  1,K,^ 


where  Fj_^  is  taken  to  be 
expression  for  F l,-  In  other  words,  for  the  Modified  Theodorsen 
Analysis , 


only  the  0  dependent  terms  of  the  classical 


K  =  1  -  1  XSY  +  AbSfl 

L  <>  r) 


and 


Mrt  =  -  -  W>^se  +  I  ALS^e. 

'^2  2 

Throughout  the  derivation  of  both  analyses,  L  was  assumed 
to  depend  only  on  the  preset  angle  of  attack  Qq  ,  for  a  given  foil  sec¬ 
tion,  and  to  be  independent  of  other  parameters  such  as  speed  S 
and  mass  unbalance  mb  .2  in  this  present  study,  these  EUialyses, 
based  on  this  assumption,  are  used  to  make  computations  of  damping 
ratio  and  frequency  of  vibration.  It  is  assumed  that  these  computa¬ 
tions  are  reasonably  accurate  so  long  as  L  does  not  vary  widely  with 
parameters  other  than  .  (See  Discussion  of  Results  and  Recom¬ 

mendations.  ) 
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APPENDIX  E 


CRITERIA  FOR  SELECTION  OF  INITIAL  CONDITIONS 
WHICH  ALLOW  ONE  MODE  TO  PREDOMINATE 
IN  THE  TOTAL  RESPONSE 


Two  modes  of  vibration  are  generally  associated  with  solu¬ 
tions  of  the  equations  of  motion  for  the  TMB  Control  Surface  Flutter 
Apparatus.  Criteria  for  the  proper  selection  of  initial  conditions 
leading  to  the  suppression  of  one  of  the  two  modes  is  now  derived. 

If  the  equations  of  motion  for  either  the  Extended  Simplified 
or  Modified  Theordorsen  Analyses, 

Bi20+  B,o(9+  BjjY  +  BjjY  =  0 


have  basic  solutions  of  the  form 


d  =  and  Y  =  Yoe^‘;  Y.e^*‘ 


where  0^,  Yg,  X  and  A*  are  generally  complex,  then  the  gen¬ 
eral  solution  of  the  differential  equations  of  motion  is 


■  If./'' 


If  the  initial  conditions  are  chosen  to  suppress  the  Xj 
mode,  solutions  of  the  equations  of  motion  representing  the  X^ 
mode  are 


0  =  + 


xTt 
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These  latter  solutions  for  0  and  Y  ,  with  one  mode  com 
posing  the  total  response,  are  also  solutions  of  a  well  known  and 
simpler  set  of  differential  equations 

+  Afl  +  iifl  =  0 


Y  +  CY  +  DY  =  0, 


To  solve  these  equations  for  0  ,  assume  0  «  c* 
Substituting  this  expression  for  0  ,  together  with  its  first  and 
second  time  derivatives,  in  the  simpler  differential  equation  in  <9 
yields  a  quadratic  equation, 


n^+  \n+H  =  0 


having  roots 


"i  = 


-  A  +  \/a  *  -  4B 
2 


and 


-A  -  Va*  -  4B 


With  expressions  for  the  roots  n,  and  n2  known,  t'\e  quadra¬ 
tic  equation,  n^  +  An  +  B  =  0  may  be  written  as 

(n  -  iij)  (n  -  Hj)  =  0 


or 


n*  —  (nj+njln+njHj  =0. 
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Comparing  this  form  of  the  equation  with  the  form  origi 
nally  used  shows  that 


\  =  -  (nj  +  02^  and  H  ain2  • 


The  general  solution  of  the  equation 


0  +  Tfl  +  llfl  =  b  -  (ttj  +  n2W  +  n|n2^  -  0 


is 


/I  r  1  ^  r  "  2  ^ 
=  r,e  +  r  e 


If  K  ,  K  ,  n,  and  nj  are  replaced  by  ^oi  •  ^02  >  ^1  and 
A*  ,  respectively,  then 

\  =  -  (A,  +  A|)  and  M  =  AjA*  . 


A  similar  procedure  may  be  followed  for  Y  .  When  that  is  done, 
it  is  seen  that 

0  =  +  <^02®^*^ 

and 

Y  _  Y  +  Y 

'  -  *  01®  ^  '  02® 


are  solutions  of  the  differential  equations 
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At  this  point,  note  that 


+  X*  {/t,  +  jwj)  +  ((ij  -  jw^)  ■=  2|l^ 


X|X*  =  (^,  +  -  jw)  =  +  (Uj  . 

I 

The  preceding  differential  equations  can  thus  be  rewritten  as 


0  -  2fi ^0  +  =  0 


Y  -  2/i,Y  ^  (/if  +  wf)Y  =  0 

The  expressions  for  0  and  Y  with  the  X2  mode  sup¬ 
pressed  must  satisfy  both  the  foregoing  simpler  differential  equa¬ 
tions  and  the  equations  of  motion  for  the  flutter  apparatus.  Hence, 
the  original  equations  together  with  the  simpler  differential  equa¬ 
tions  may  be  used  to  eliminate  the  second  derivatives  and  to  yield 
a  first-order  differential  equation  which  must  be  satisifed  by  0,  Y 
and  ft ,  \  The  elimination  of  the  second  derivative  gives  the  foL 
lowing  set  of  equations  which  are  the  resulting  restriction  on  the 
initial  conditions: 

(B,i  +2,i,Bi2)0  +  [Bio-(,if  +  wf)B,2]^  + 

+  (B2i  +  2/iiB„)Y-(#if  +  oif)BjjY  =  0 


(Bn  +2/iiRj2)fl  +  [Rn,  -  (/if  +c.if)R,2l0  + 


+  (R,,+2/i,R,2)Y-[R,„-(/if  +  a,f)R,j]Y  =  0 

These  equations  enable  us  to  express  any  two  of  the  quanti¬ 
ties  ,  f)  .  Y  ,  or  Y  in  terms  of  the  other  two.  Because 

these  equations  hold  at  all  values  of  time,  they  hold  at  t  =  0  . 

Thus  these  two  equations  constitute  the  me2Uis  whereby  initial  condi¬ 
tions  may  be  chosen  to  suppress  one  mode  in  the  total  response. 

That  is.Ylf^pand  '?]t  =  o  which  may  be  arbitrarily  chosen,  and  the  cor¬ 
responding  values  Yjj^gand  01^ ^ ^computed  from  the  above  equations, 
are  the  initial  conditions  required  in  solving  the  original  differential 
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equations  (practically,  by  means  of  an  analog  or  digital  computer). 
There  is  then  no  contribution  of  the  Xj  mode  to  0  (t)  and  V  (t). 


APPENDIX  F 

THE  EFFECT  OF  INITIAL  CONDITIONS  ON 
THE  VIBRATION  OF  A  SHIP'S  RUDDER 


Solutions  Y  (t)  and  ^(0  of  the  equations  of  motion  for 
the  flutter  apparatus  represent  the  superposition  of  two  modes  of 
vibration  when  arbitrary  initial  conditions  are  assumed.  The  solu¬ 
tion  can,  how4v§r,  be  restricted  to  represent  a  single  mode  of 
vibration;  i.  e.  ,  one  of  the  two  modes  cam  be  suppressed  if  the  initial 
conditions  are  properly  chosen.  This  may  be  seen  from  the  analog 
solution  shown  in  Figure  3.  Similar  solutions  can  be  obtained  by 
the  digital  method  of  computation.  Criteria  for  the  proper  selection 
of  these  initial  conditions  are  established  in  Appendix  E. 

Similarly,  the  vibratory  response  of  any  point  on  a  rudder 
or  movable  control  surface  of  a  ship  will  also  depend  upon  the  initial 
conditions  of  the  rudder'^*’’  in  translation  and  rotation.  In  actual 
practice  the  initial  conditions  of  a  rudder  vibration  may  be  estab¬ 
lished  by  considering  an  impulsive  force  of  duration  tQ  applied  at 
a  given  point  on  the  rudder  at  time  t  =  0.  The  response  of  the  sys¬ 
tem  immediately  after  the  imposition  of  this  force  (after  t  =  to)  is 
equivalent  to  ftie  solution  of  the  homogeneous  differential  equations 
of  motion  for  the  rudder  system  with  initial  conditions  the  same  as 
the  condition  of  the  forced  system  at  t  =  tQ- 

To  obtain  the  actual  initial  conditions  at  time  t  =  tQ  for  the 
homogeneous  equations,  it  is  necessary  to  solve  the  inhomogeneous 
equations  of  motion  which  include  the  impulsive  force.  Solutions 
will  express  the  motions  of  the  rudder  in  translation  and  rotation  as 


’’‘This  is  easily  done  on  the  analog  computer.  Note  that  all  the  co¬ 
efficients  are  considered  as  known,  /j,  and  cuj  ,  having  been  pre¬ 
determined  on  the  analog  or  digital  computer  as  previously  explained. 

**The  term  rudder  will  hereafter  be  used  in  a  broad  sense  to  include 
rudders  and  all  rudderlike  appendages  of  ships,  such  as  submarine 
diving  platnes. 
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functions  of  time.  The  numerical  values  of  these  quantities  and  their 
time  derivatives,  evaluated  at  t  »  are  the  initial  conditions 

of  the  vibrating  rudder  system. 

It  may  be  inferred  from  Table  1  of  Reference  9  that  the 
three  modes  of  vibration  associated  with  athwartship  motions  for  an 
actual  rudder  system  will  in  general  exist  simultanequsly,  unless 
the  initial  conditions  are  specially  chosen.  Contributions  to  the  dis¬ 
placement  of  any  point  on  the  rudder  for  each  coupled  mode  are  due 
to  the  motions  in  the  translational  and  the  two  rotational  degrees  of 
f  reedom. 


The  influence  of  a  given  set  of  initial  conditions  on  the  vi¬ 
bratory  response  at  different  points  on  the  rudder,  measured  simul¬ 
taneously,  vaII  generally  vary  because  the  initial  conditions  will 
make  different  modes  (and  their  associated  frequencies)  predominate 
in  the  vibratory  response  at  different  points  on  the  rudder.  In  other 
words,  different  predominant  frequencies  may  be  measured  at  dif¬ 
ferent  points  on  a  rudder  at  the  same  time.  In  addition,  the  relative 
vibratory  response  at  those  same  points  will  vary  with  chauiges  in 
the  initial  conditions,  or  its  equivalent,  the  point  of  application  of 
the  impulsive  force. 

Finally,  at  a  given  point  on  a  rudder  the  frequency*®*  of  a 
given  response  curve  may  vary  if  an  initially  predominant  mode 
decays  more  rapidly  than  an  initially  imperceptible  mode. 

The  previous  discussion  of  initial  conditions  is  substanti¬ 
ated  by  the  analysis  of  the  rudder  vibration  records  for  ALBACORE 
presented  in  Reference  10. 

APPENDIX  G 
ANAUXIY  BETWEEN 

TMB  CONTROL  SURFACE  FLUTTER  APPARATUS 
AND  THE  RUDDER  CONFIGURATION 
OF  USS  FORREST  SHERMAN  (DD  931) 

The  following  comments,  comparing  experimental  condi¬ 
tions  for  the  flutter  apparatus  with  trial  conditions^  of  FORREST 
SHERMAN,  ^  (on  which  apparent  flutter  phenomena  of  the  rudders 
was  observed)  were  made  by  D.  A.  Jewell  in  Reference  7.  These 
comments  are  considered  to  be  of  great  interest  in  connection  with 
the  present  study  and,  because  of  the  concise  wording  of  the  original, 
are  quoted  here. 

*^more  exactly,  the  interval  between  axis  crossings. 
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"in  many  ways,  the  experimental  conditions''  were 
similar  to  those  on  the  DD  931.  The  two  degrees  of  freedom, 
the  natural  frequencies,  the  general  rudder  configuration,  and 
location  of  the  rotational  axis  with  respect  to  the  rudder  were 
the  same  in  both  cases.  Two  or  three  experimental  parametric 
values,  however,  differ  significantly  from  DD  931  values.  In 
dimensionless  form,  these  parameters  are: 


The  mass -density  ratio  - 


/'  ^ 


in 


The  distance  of  the  center  of 
mass  aft  of  the  rotational  axis  -  . 


The  frequency 


k 


o)  b 

V 


where  m 

P 

b 

I 

<0 


is  the  translational  mass, 
is  the  fluid  density, 
is  the  rudder  semichord, 
is  the  rudder  span  length,  and 
is  the  frequency  at  speed  V. 


VALUES  OF 

THREE  DIMENSIONLESS  PARAMETERS 


Parameter 


Experimental 


DD  931  Class  Ships 


P 


=  3.  0 


36 


X 


a 


0.  23  to  0.  31 


0.  20 


k 


=  1.0 


2.  3 


The  different  values  of  fi  and  x„  for  the  ship 
strongly  indicate  that  the  ship's  critical  flutter  speed  is  much 
higher  than  the  experimental  speeds.  For  instance,  at  ' 

*a  =  0-31  ,  where  p  was  increased  from  3.  24  to  3.  56 
(configuration  B  to  D),  the  critical  flutter  speed  increased  from 

^  Configurations  B  and  D. 
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approximately  11.  5  to  13.  5  knots.  The  effect  of  the  different 
values  of  L  as  an  independent  parameter  cannot  yet  be 
evaluated. 

The  usefulness  of  the  Modified  Theodorsen  Analysis  in 
predicting  the  stability  of  the  flutter  apparatus  can,  therefore, 
reasonably  be  regarded  as  an  indication  of  the  usefulness  of  that 
analysis  in  predicting  ship  control- surface  flutter,  after  the  refine¬ 
ments  suggested  in  the  Recommendations  have  been  made. 


APPENDIX  H 

THE  CLASSICAL  THEODORSEN  EQUATIONS 
AND  THEIR 

QUESTIONABLE  APPLICABILITY  TO  RUDDERS 


The  classical  Theodorsen  equations  are  now  written  for 
the  case  of  a  thin  flat  foil  of  width  2b  and  infinite  length,  inclined 
at  an  angle  a  to  a  uniform  stream  of  velocity  v;  see  Figure  88. 


Let  hg  be  the  displacement  of  the  foil  centerline  trans¬ 
verse  to  V  and  let  the  foil  be  stationary  or  execute  small  vibrations 
at  circular  frequency  w  .  Then  the  reactions  on  the  foil  are  equiv 
alent  to  three  distributed  forces  and  a  distributed  moment  whose 
respective  magnitudes  per  unit  length  are:’'* 

*  This  formulation  can  be  obtained  from  equations  XVIII  and  XX  of 
Reference  5  by  omitting  all  T  terms  (no  aileron  considered  here), 
inserting  C(k)  =  F  +  iG  and  making  everything  real  by  writing 
a  «  ^  »  and  »<i  =  ^  ;  and  similarly  for  h;  let  h  =  -hQ. 


Also  write  P  =  -(Pj  +  Pz  +  P3)  and  choose  a  =  0  so  that  the  axis  of 
rotation  is  the  midchord  line.  Then  note  that  if  Pj,  P2,  and  P3  are 
assumed  to  act  along  the  lines  specified  above,  they  provide,  in  the 

total  moment  about  the  midchord  line,  all  terms  except  that  called 
M  above.  That  is,  with  the  above  substitutions,  equation  XX  gives 


Mfl  =  ~  r  Trph*a  -  ffpb*vd  +  ffpvb^(F  +  iG)(v  a  -  hg  +  bd) 

8  2  2 

=  -  wpb*  »fpb*  vd  +  npvh^  F(va  -  bg  +  -^  ba)  +  np\b^  (vd  -  bg  +  -^  b&’)  • 

8  2  2  w  2 


In  the  notation  used  here,  =  M  +  -^  .  Other  lines  of 

action  could  be  assumed  for  Pj,  P2,  and  P3;  then  the  expression  for 
M  would  have  to  be  corrected  so  as  to  keep  the  value  of  the 

same. 
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Figure  88 

Forces  and  Moment  on  Foil  Moving  in  a  Uniform  Stream 
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=  2»Tpbv*^ - +  irpb*v  (K  +  Y  G)a  -  2ffpb*  ^  G(hQ  -  j  a  ) 


1%  =  ffpb 


Pj  =  -  npb^hg 

M  =  —  i  ffpb*  a 
8  ^ 

Here 


acting  along  the  forward  quarterline  A, 
acting  along  the  rear  quarterline  B,  and 
acting  along  the  midchord  line,  and 


V 

bo) 


4I-  ;  T  =  period  = -^ 
>r(2b)  (o 


F  and  G  are  functions  of  7  (see  Reference  5,  Figure  4), 
with  F  positive  and  G  negative.  F  rises  from  F  =  ^  at  7  =  0  to 
1(?)  '■  as  -  oc  .  G  decreases  from  G  =  0at  7  =  0toa  minimum 
of  -0.  20  at  about  (?)  ^  =  5  and  then  slowly  increases  to  0  (?)  as 


can  be  visualized  as  the  distance  that  the  stream  ad¬ 
vances  during  one  period  of  vibration  divided  by  n  times  the  foil 
width. 


In  these  formulas,  a  represents  rotation  of  the  foil  about 
the  midchord  line.  It  maybe  preferable  to  let  the  rotation  a  occur 
about  a  parallel  axis  e  away,  e  being  positive  ahead  of  midchord. 
Let  bg  denote  the  displacement  of  the  foil  at  this  new  axis;  then 
I’d  =  *'e  “  •  Substitution  for  in  the  formula  for  Pj 

changes  the  first  two  terms  to  Pj  where 

Pj  =  2ffpbv*  F  +  ffpb^v  ^  ^  T  ■ 

The  individual  magnitudes  of  these  terms  are  thereby  changed,  but 
their  sum  remains  the  same.  It  is  simplest,  however,  to  let  hQ 
stand  in  the  third  term  of  Pj  and  in  P3  and  M.  The  total  turning 
moment  Mg  about  the  new  axis  of  rotation  is  given  by  the  formula 

M.  -  P.(|-«)-P.(|  +•)  -Pi'+M 

♦These  (? )  are  usedbecause  the  graphs  go  only  to  -^  =  40;  values  at 
7  "““have  not  been  verified;  or  graph  is  on  too  small  a  scale. 
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The  forces  that  give  rise  to  M  are  distinct  from  those  in¬ 
cluded  in  Pj,  P2.  or  P3;  they  add  up  to  a  zero  resultant  but  provide 
a  turning  moment  M,  of  the  magnitude  stated,  about  any  axis, 

DISCUSSION 

The  validity  of  Theodorsen's  formulas  for  rudders  is  dis¬ 
cussed  presently.  If  his  formulas  are  used,  it  appears  probable 
that  the  constant  ^  will  normally  be  less  than  unity  for  ships.  When 

■^  <  I  ,  ^  K  <  O..Ti> ,  —  (>  <  0.0^  , 

Lift  Force  If  fi  -  0  and  =  0,  only  the  first  term 
of  Pj  remains  as  the  familiar  lift  force  per  unit  length  for  steady 
motion  acting  along  the  forward  quarter- line.  Here  hg  represents 
a  steady  translation  and  merely  alters  the  steady  angle  of  attack.  In 
this  case  w  =  0,  hence 

IT  ^  ^  and  (presumably)  F  =  1. 

Virtual  Mass  P3  and  M  may  be  regarded  as  reac¬ 
tions  to  a  distributed  virtual  mass  m^  along  the  midchord  line  and 
a  distributed  virtual  moment  of  inertia  along  this  line.  The  coef¬ 
ficient  of  (hg  -  -ja)  in  Pj  will  be  denoted  by  Mp  .  The  magni¬ 
tudes  are  (per  unit  length): 

Mp  =  nph^  ;  =  -  2*rpb*  5  V  *  )  . 

•  •  L 

Note  that  hg  -  f  o  is  the  acceleration  of  the  rear  quarter  line  and 
that  G  S  0,  so  that  M„  2  0.  Also,  —  =  .  The  quantities  \1„  and  L 

^  ci>  k  f'  r 

have  the  usual  values  for  a  uniform  plane  strip  of  infinite  length 
moving  transversely. 

Significant  Terms  in  Pi  The  quantity 

L  <  2,  Mp/Mp  <  0.6.  As  <u-»0  and -•«,Mp -•'»(?)  but  less  rapidly  than  l/cu; 
however,  the  contribution  of  the  term  to  Pj  for  given  amplitudes 
of  hg  and  a  is  numerically  proportional  to  bv  w  G,  and  so  -*  0 

as  fo  ^  0  (since  =  -<u'^hg,a  =  -w^a  ).  This  part  of.Pj  is  not  enor¬ 
mously  large,  as  shown  by  a  comparison  with  the  first  two  F  terms 
in  Pj.  The  terms  in  h'o  and  Iiq  differ  in  phase,  but  may  be  compared 
as  to  magnitude  by  introducing  the  amplitudes  Kq  amjr  amplitude  of 
Ho,  etc.  ;  then  Hq  amp  ‘^Hpamp,  Similar  comparisons  may 

be  made  for  a  and  a  .  Then  it  is  seen  that  (since  J-  = 

y  k  DO)  y 
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Iiq  ratio 


'Inp  {bv/(/j)  (  -  (1)  li,| 
2>rpbvl''lio 


£ 


2ffp  (bv/w)  ( -(i)  5  (bw^a,„p) 
2rrpbv2Ka.„p 


1 

2 


Here 

—  <0.3!>(?);  -  1  <  0.20(?). 

K  2  \  F  / 

•  •  . 

Hence  for  a  rough  approximation,  the  G(ho  -■|a)tc>.'m  in  P  perhaps  may 
be  ignored  in  comparison  with  the  first  term  of  P^. 


Damping  Terms  Two  methods  for  discussing  the 
foil  damping  are  now  presented.  In  the  first  method  the  forces  and 
moments  on  the  foil  associated  with  the  a  terms  are  considered, 
whereas  in  the  second  method  the  energy  of  the  foil  associated  with 
the  Kq  and  a  terms  is  treated. 


Force- Moment  Method  The  d  terms  in  P2 
and  Pj  add  to  a  force  P^  of  magnitude 

=  nph  (1  +  F  +  G)  d  . 

Also,  these  two  terms  in  P2  and  P^,  one  acting  at  -  and  the 
other  at  +  j  ,  give  rise  to  a  turning  moment  M^  about  midchord 
of  magnitude 


=  -i  ffpb®  V  ( - 1  +  F  +  -^  G)d  . 

M^/d  <  0  since  F  <  1  and  G  <  0.  Our  equations  are  valid  only 

for  harmonic  motion,  but  we  surmise  that  whatever  physical  process 
gives  rise  to  a  ratio  ^^tj/d  that  is  negative  in  all  cases  of  har¬ 

monic  motion  will  tend  to  damp  out  any  motion,  i.  e.  ,  have  a  positive 
damping  effect.  The  effect  of  P^j  ,  however,  is  not  easily  seen. 

A  rough  estimate  of  the  quantitative  signijficance  of  P^; 
and  M(i  may  perhaps  be  made  by  comparing  these  quantities,  re¬ 
spectively,  with  the  a  term  in  P.  aund  with  the  turning  moment 
about  midchord  due  to  this  term.  Because  of  a  90  deg  difference  in 
phase,  the  comparison  is  made  most  easily  in  terms  of  amplitudes 
by  recalling  that  «  a.mp  ~  amp’  algebraic  signs  are 

then  of  minor  interest;  hence,  for  convenience  -M^  .  will  be  used 

instead  of  +M  ^  Denoting  the  amplitude  ratios  for  P  and  M  by 
rp  and  rj^  respectively: 
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= 


ffpb*v(l  +  K  +  "J"  G)wn,mp 
2»rpbv*Ka.„p 

-iffpb»v(l-  b’-IOwa.^p 


I’M  = 


ffpb^v^  Fa 


amp 


c 

F 


G 

V 


since  b^  ^  k.  Some  values  calculated  from  these  formulas,  using 
Theodorsen's  Fig.  4  for  F  and  G,  are: 


s 

0 

0.  25 

B 

1.  0 

2.  0 

4.  0 

10.  0 

20.  0 

30.  0 

40.  0 

00 

F 

5 

m 

0.  52 

0.  54 

0.  70 

0.  84 

0.  91 

0.  94 

0.  96 

1.  00 

G 

1 

-0.  06 

-0.  15 

-0.  18 

-0.  16 

-0.  13 

-0.  10 

0.  09 

0 

1 

5.  9 

2.  80 

1.  24 

0.  42 

0.  05 

-0.  08 

-0.  09 

-0.  07 

-0.  07 

0 

00 

1.  98 

1.  04 

0.  61 

0.  42 

0.  31 

0.  20 

0.  15 

0.  11 

0.  09 

0 

From  these  values  of  rp  and  rjy^,  it  appears  that  when  1/k  <0.  25, 

at  least,  the  d  terms  in  Pj  and  may  be  expected  to  have  large 
effects  in  comparison  with  the  a  term  in  P^.  Only  when  1/k  >  5 

do  the  terms  become  relatively  rather  small. 

Energy  Method  The  terms  containing  hQ  or 
d  in  Pj  and  P2  have  the  form  of  damping  terms.  It  is  worthwhile 
to  ascertain  the  sign  of  their  net  effect  upon  the  energy  of  the  vibrat¬ 
ing  foil. 


Let  dW/dt  stand  for  the  inflow  of  energy  into  the  foil  due 
to  the  components  of  force  represented  by  the  Kq  euid  d  terms.  If 
dW/dt  <  0,  the  net  energy  flow  is  outward  as  in  ordinary  (positive) 

damping. 


The  value  of  dW/dt  can  be  found  by  multiplying  each  com¬ 
ponent  of  force  by  the  velocity  of  its  point  of  application,  and  adding 
these  products.  The  velocity  is  hg  +  bd/2  for  any  term  of  Pj,  but 
Kq  -  bd/2  for  P2»  which  acts  along  the  rear  quarter- chord  line. 
Hence  (in  the  case  of  harmonic  motion): 
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(It 


-  ■  -  - W  *  *  I ■')  -  l'‘(s  -  ?'‘)j 


or 


(J« 

(It 


I  -  F  - 


The  factor  enclosed  in  brackets  [  ]  in  the  last  equation 
will  be  denoted  hereafter  by  the  symbol  [  1  .  The  Iiq  and  d  ^ 

terms  are  positive,  since  F  <  1  and  G  <  0,  but  the  sign  of  the 

h'o  d  term  is  not  fixed.  Hence,  a  numerical  investigation  is  nec¬ 

essary. 


Write  ^  =  flho  where  6  may  have  any  real  value.  Then 
after  dividing  by  h* 

4  [  1  =  F  -  (I  +  4  +  (1  -  F  -  I-  G)e\ 

n  Q  k  k 

If  <?  -•  +  oo  ,  the  0^  term  predominates  and  [  ]  must  be  posi¬ 
tive,  since  F  <  1  and  G  <  0.  Between  these  extremes,  [  l/h^, 

regarded  as  a  function  of  &  ,  must  have  a  minimum  value  and 
[  ]  can  become  negative  only  if  this  minimum  is  itself  negative. 

• 

To  locate  the  minimum,  set  the  derivative  of  [  ]/h* 
with  respect  to  0  equal  to  zero  and  solve  for  d  ,  obtaining 


e  = 


1  +  1 


G 


2(1  -F  -  -|G) 


Substitution  of  this  value  of  9  gives  as  a  minimum  value 


F  - 


(1  +-|G)^ 
4(1  -  F  -  f  G) 


Some  values  calculated  with  the  use  of  Theodorsen’s  Figure  4  are; 
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It  may  be  concluded  that,  at  least  for  0  =  1/k  <  21, 

I  I  '  0  and  hence  dW/dt  ^  0,  so  that  the  net  effect  of  the  hg  and  « 

terms  in  Pj  and  P2  is  to  withdraw  energy  from  the  foil.  Thus,  if 
steady  flutter  occurs  with  1/k  in  the  range  specified,  other  terms 
must  provide  an  equal  inflow  of  energy.  At  speeds  below  the  lowest 
speed  at  which  flutter  is  possible,  it  may  be  surmised  that  the  physi¬ 
cal  processes  responsible  for  the  lig  and  (i  terms  in  the  formulas 
for  the  harmonic  case  will  give  rise  to  positive  damping. 

It  would  seem  that,  if  the  h^  term  in  Pj  is  retained,  as  is 
usually  done,  then  for  consistency  the  'i  terms  should  be  retained 
also.  In  steady  motion  the  hg/v  term  results  from  a  simple  cause  - 
namely,  a  change  in  the  angle  of  attack  -  whereas  the  origin  of  the 
fi  terms  is  not  obvious,  but  this  difference  scarcely  seems  to 
justify  retention  of  the  hg  term  alone. 

Simplified  Approximations  The  following  approxi¬ 
mations  may  be  useful  for  some  purposes.  The  steady-motion  lift 
force  described  on  page  836  is  frequently  used  as  an  approximation 
even  when  n  is  allowed  to  vary.  It  is  then  necessary  to  consider 
whether  Hg  should  be  replaced  by  the  displacement  of  some  point 
other  than  the  forward  quarter-line,  and  what  assumption  should  be 
made  concerning  the  effective  line  of  action  of  the  lift  .  If  this  line 
of  action  is  distant  I.  from  the  axis  of  rotation,  the  turning  moment 
about  this  axis  will  equal  I.  times  the  lift  force. 

For  a  closer  approximation  the  d  terms  may  be  retained. 
So  long  as  1/k  does  not  exceed  1  or  2,  however,  the  error  may  not 
be  excessive  if  all  G  terms  are  omitted  and  F  is  replaced  by  1/2. 
This  was  done  in  Reference  2  in  obtaining  the  "Modified  Theodorsen 
Analysis.  "  To  obtain  the  results  in  Reference  2,  we  assume  that  the 
originally  specified  lines  of  action  for  Pj  and  P_  are  retained.  If 
F  =  1/2  and  G  =  0,  the  results  for  PJ  and  P2  wifti  the  axis  of  rotation 
distant  "e"  ahead  of  midchord  are: 

-i'j.  .pb>v 

Pj  =  npb^v  d  . 

Therefore,  the  total  lift  force  per  unit  length  on  the  foil  is: 

+  -^)d. 

v/  \2  b/ 


Pa  = 


Pj  +  P2  =  ’rpbv 
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We  do  not  try  to  specify  a  line  of  action  for  P-^  as  a  whole. 
To  calculate  the  total  moment  Ma  about  the  axis  (exclusive  of  in¬ 
ertial  moment  M),  it  is  necessary  to  split  again  into  its  two  parts, 
Pj  acting  at  (b/2  -  e)  ahead  of  the  axis  and  P2  acting  at  {-b/2  -e)  behind 
the  axis  (e  is  negative  if  behind  midchord).  Then 


Mg  =  p;  (b/2  -  e)  +  Pa(-b/2  -  e)  =  P;  (b/2  -  e)  -  P2(b/2  +  e) 

■»  ffpbv*(a  -  hg/v)(b/2  -  e)  +  wpb*vd(l/2  +  e/b)(b/2  -  c)  -  ffpb*vd(b/2  +  e) 
=  -ffpbv*(b/2  -  e)  (a  -  b^/v)  +  ffpb*v(b/4  +  c/2  -  e/2  -  e*/b  -  b/2  -  e)d 
=  irpbv*(b/2  -  e)  (a  -  hg/vl  -  ffpbv(b/2  +  e)*d  . 


In  the  special  case  e  =  b/2: 


“  ffpbv*(a  -  hj/v)  +  2jTpb*vd 
Mg  =  -  wpb*  vd  . 

These  values  of  Pg,  and  Ma  agree  with  those  of  Fj^  and 
as  given  in  the  middle  of  page  37  of  Reference  2,  where  S  =  v, 

Y  -  h_,  0  =  a  ,  arid  A  =  Znphl.  The  length  of  the  foil  I  occurs  in 

A  because  F]_,  and  represent  total  reactions  on  the  foil,  whereas 
Pg^  and  Ma  stand  for  magnitudes  per  unit  length. 

No  specific  line  of  action  for  Pg  or  F^  is  assumed  here, 
since  the  total  turning  moment  about  the  a  or  9  axis  is  explicitly 
stated  in  the  second  equation.  That  is,  the  same  distributed  forces 
whose  resultant  is  Pg^  or  F^  give  rise  to  the  moment  M^  or  M^. 
The  usual  device  of  replacing  the  second  equation  by  Ma  =  LPg 
or  Mj  =  L  F^  is  not  available  because  M  a  or  M^  is  not  propor¬ 
tional  to  the  whole  of  Pg  or  F^.  To  facilitate  thinking,  however,  it 
may  be  desirable  to  assume  that  Pg  or  Fj^  acts  along  a  line 
through  a  or  9  axis. 
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DAMPING  COMPUTATIONS  FOR  ALBACORE  RUDDER 
AND  TMB  FLUTTER  APPARATUS 

A  meaningful  interpretation  of  rp  and  r j.  is  strictly  valid 
for  the  foil  described  on  page  83^  and  approximately  valid  for  an  air¬ 
plane  wing;  however,  for  a  structurally  damped  rudder  with  finite 
aspect  ratio  and  thickness,  attached  at  the  quarter  chord(e  ^  0) 
through  a  flexible  rudder  shaft  to  a  flexible  hull,  the  actual  values  of 
rp  and  rj^  would  undoubtedly  be  different  from  the  theoretical  values 
found  from  the  equations  derived  on  page  838.  Moreover,  the  motions 
of  the  top  of  the  rudder  stock,  at  the  point  of  attachment  of  the  rudder 
stock  to  the  hull,  certainly  affect  the  damping  of  the  rudder.  Finally, 
the  athwartship  vibratory  displacements  of  the  rudder  consisting  of 
motions  in  three  degrees  of  freedom,  one  in  translation  and  two  in 
rotation,  further  invalidate  the  use  of  rp  and  rj^  in  determining  the 
significance  of  the  P^  and  M^  terms. 

For  the  TMB  Flutter  Apparatus,  the  relatively  rigid  rudder 
stock  (or  shaft)  is  attached  at  its  upper  end  to  a  rigid  structure  and 
the  hydrofoil  motions  are  restricted  to  two  degrees  of  freedom. 
However,  again  e  0,  the  foil  is  structurally  damped  and  has  a  fi¬ 
nite  aspect  ratio  and  thickness. 

Nevertheless,  it  is  instructive  to  show  by  example  how, 
based  on  Theodorsen's  equations,  a  value  of  r  and  rj^  can  be  deter¬ 
mined  for  an  actual  or  rudder  model  (i.  e.  ,  rudder  plus  stock  at¬ 
tached  to  a  rigid  hull)  when  the  effects  of  finite  aspect  ratio  and 
thickness  are  ignored  and  when  the  motions  of  the  rudder  in  onlytwo 
degrees  of  freedom  are  treated.  The  motions  in  these  degrees  are 
considered  to  make  the  most  significant  contributions  to  the  vibra¬ 
tory  displacement.  Although  the  results  here  for  rp  and  rj^  hold 
only  for  e  =  0,  the  procedure  for  deriving  the  equations  leading  to 
a  single  calculation  for  rpand  rjy^  points  the  way,  for  the  interested 
reader,  towards  deriving  just  as  easily,  similar  equations  for  "e" 
corresponding  to  a  rudder  turning  about  a  rudder  stock  attached  at 
the  quarter-chord,  i.  e.  ,  e  =  • 

2 

For  the  rudder  of  ALBACORE,  b  may  be  estimated  as  3  ft. 
and  the  lowest  vibration  frequency  for  \ya  vibrations  as  780  cpm’r 
Thus  rj  =  2  rr  X  780/60  =  82  rad/ sec  and,  if  v^^^is  the  speed  in  knots. 


*As  shown  in  this  report,  the  frequency  does  not  vary  much  with 
speed. 
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1/k  =  v/b^j  =  1.69v,,„  (3x82)  - 


=  10  20  30 

1/k  =  0.060  0,138  0.207  , 

At  such  speeds  the  chart  on  p.  838  indicates  that  rp  >  6,  rj^  ^  2, 

Such  values  of  r  and  rj^  discourage  omission  of  the  a  terms  in 
Pj  and  P2 

For  the  TMB  Flutter  Apparatus  described  in  Reference  2, 
b  =  9  in.  and  the  frequency  is  4  cps,  hence 

l/k  =  11.69vk„/0.78/(tx2>7)l  =  Vk„/ll  . 

Thus  at  22  knots,  l/k  =  2  and  rp  =  0,  42,  rj^  =  0.  42.  Conceivably 
such  values  of  rp  and  r  justify  the  suspicion  that  failure  to  obtain 
flutter  in  the  observations  might  have  been  due  largely  to  excessive 
hydrodynamic  damping. 

Perhaps  for  airplane  wings,  low  fj  and  high  v  make  k 
and,  hence,  also  r  and  r.^  much  less  than  1,  so  that  neglect  of  the 
fi  terms  is  justified, 

APPUCATION  OF 

THE  CLASSICAL  THEODORSEN  EQUATIONS 
TO  RUDDERS 

Theodorsen’s  results  refer  to  the  case  of  a  uniform  flat 
foil  of  infinite  length.  In  contrast,  the  aspect  ratio  of  ships'  rudders 
(height  over  width)  is  rather  small  (always  between  1  amd  3).  The 
forces  in  this  case  are  likely  to  be  decidedly  different  from  those 
obtained  by  direct  use  of  the  classical  Theodorsen  expressions. 
Recent  work  described  in  various  reports  does,  in  fact,  indicate 
that  this  difference  is  probably  large.  (See  discussion  by  H,  N. 
Abramson  in  Reference  4.  ) 

When  the  foil  is  locked  in  rotation  ( d  =  a  =  0)  at  an  angle 
of  attack  a  =  ,  then  for  a  constant  displacement  hg  (ho  ^*0  ~ 

the  steady  lift  force  (c(k)  =  1)  given  by  Theodorsen's^classical  equa¬ 
tion  is  (for  a  foil  of  length  I  ): 

Pj  =  Pjl  =  2  77’p/bv^aQ  • 

Whereas  the  Extended  Simplified  or  Modified  Theodorsen  Analyses 
give  (In  Appendix  D  let  Fj_^,  S,  and  d  be  denoted  by  Pj ,  v,  and  a 
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respectively.  ) 


P',  =  . 

For  the  TMB  Flutter  Apparatus,  it  was  found  that  the  lift 
force  predicted  by  the  classical  Theodorsen  analysis  is  approxi¬ 
mately  three  times  as  large  as  the  lift  force  found  by  (experiment;^* 
i.  e.  ,  the  lift  constant  Inp  I  h  =  .'t  A  ,  where  the  experimental 
value  for  A  is  0.034  lb  -  sec^  (See  Table  2.  )* 

in.  2 

Actually,  not  much  is  known  so  far  about  the  magnitude  of 
the  force  and  moment  reactions  on  a  ship's  rudder,  except,  of  course, 
in  the  case  of  a  steady  rudder  angle.  In  other  cases,  measurement 
of  the  reactions  is  no  doubt  difficult.  Attention  should  be  directed  at 
present  toward  obtaining  actual  measurements  of  the  reactions  on 
vibrating  rudders  or  rudder  models;  such  experimental  data  may  lead 
to  the  establishment  of  a  reliable  method  of  predicting  both  steady 
and  unsteady  lift  force  and  moment  vs.  speed  relations. 

Moreover,  as  pointed  out  in  the  Discussion  of  Results,  the 
inclusion  or  omission  of  structural  damping  can  make  a  relatively 
large  difference  in  the  predicted  critical  flutter  speed. 

Free  surface  waves  and  cavitation  affect  the  boundary  con¬ 
ditions,  thereby  causing  additional  differences  between  hydroelastic 
euid  aeroelastic  phenomena. 

Further  reservations  as  to  the  applicability  of  Theodorsen's 
equations  to  ship  and/or  model  rudders  are  based  on  analysis  and 
evaluations  given  in  Discussion  of  Results. 


Note  that  2rrplb  =  6.  28  (9.  34  x  10”^  lb-sec^)  (18  in.  )  (9  in.  )  = 

in^ 

0.  094  lb  -  sec^ 

in  2 
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Theodorsen  Analysis  gives  a  consistently  better  agreement  Theodorsen  Analysis  gives  a  consistently  better  agreement 

with  experimental  data  than  does  the  Extended  Simplified  with  experimental  data  than  does  the  Extended  SimpliHed 

Analysis;  moreover,  the  results  of  the  former  analysis  are  in  Analysis;  moreover,  the  results  of  the  former  analysis  are  in 

good  agreement  with  available  experimental  data,  whereas  the  good  agreement  with  available  experimental  data,  whereas  the 

results  of  the  latter  analysis  are  not.  results  of  the  latter  analysis  are  not. 


